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Abstract 
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M. Zukowski and D. Kaszlikowski, Vienna Circle Yearbook vol. 7, editors D. Green- 
berger, W. L. Reiter, A. Zeilinger, Kluwer Academic Publishers, Dortrecht (1999) 



[0 M. Zukowski, D. Kaszlikowski, A. Baturo and J. -A. Larsson, |quant-ph/991005^ 



D. Kaszlikowski, P. Gnacinski, M. Zukowski, W. Miklaszewski and A. Zeilinger, [quant- 
ph/0005028| 



The first two chapters have an introductory character. 

In Chapter 3 it is shown that the possibility of distinguishing between single- and 
two- photon detection events, usually not met in the actual experiments, is not a necessary 
requirement for the proof that the experiments of Alley and Shih ||Shih88|| and Ou and Mandel 



L)u88|| are, modulo a fair sampling assumption, valid tests of local realism. It is also shown 



that some other interesting phenomena (involving bosonic-type particle indistinguishability) 
can be observed during such tests. 

Next in Chapter 4 it is shown again p| that the possibility of distinguishing between single 
and two photon detection events is not a necessary requirement for the proof that recent 
operational realisation of entanglement swapping cannot find a local realistic description. 
A simple modification of the experiment is proposed, which gives a richer set of interesting 
phenomena. 



^The author's papers will be cited here and in the main text using the number from the list in this page 
in square brackets, whereas other papers will be quoted by the name of the first author and the year of 
publication. 
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In Chapter 5 a sequence of Bell inequalities for M particle systems, which involve three 
settings of each of the local measuring apparatuses, is derived 0]. For Greenberger-Horne- 
Zeilinger states, quantum mechanics violates these inequalities by factors exponentially grow- 
ing with M. The threshold visibilities of the multiparticle sinusoidal interference fringes, for 
which local realistic theories are ruled out by these inequalities, decrease as (2/3)*^. 

In Chapter 6 the Bell theorem for a pair of two two-state systems (qubits) in a singlet 
state for the entire range of measurement settings is presented |Q. 

Chapter 7 is devoted to derivation of a series of Greenberger-Horne-Zeilinger paradoxes 
for M quA^its (particles described by an N dimensional Hilbert space) that are fed into M 
unbiased 2A^-port spatially separated beam splitters If. 

In Part III a novel approach to the Bell theorem, via numerical linear optimisation, is 
presented J^, ||]. 

The two-qubit correlation obtained from the quantum state used in the Bell inequality is 
sinusoidal, but the standard Bell inequality only uses two pairs of settings and not the whole 
sinusoidal curve. The highest to-date visibility of an explicit model reproducing sinusoidal 
fringes is -. We conjecture from a numerical approach |^ presented in Chapter 8 that 
the highest possible visibility for a local hidden variable model reproducing the sinusoidal 
character of the quantum prediction for the two-qubit Bell-type interference phenomena is 
In addition, the approach can be applied directly to experimental data. 

In Chapter 9 the approach presented in Chapter 8 is applied to three qubits in a maximally 
entangled Greenberger-Horne-Zeilinger state. For the first time the necessary and sufficient 
conditions for violation of local realism for the case in which each observer can choose from 
up to 5 settings of the measuring apparatus are shown. 

In Chapter 10 using the modified approach developed in Chapter 8 it is shown that vio- 
lations of local realism are stronger for two maximally entangled quA^its, than for two qubits 
[§]. The magnitude of violation increases with A^. It is objectively defined by the required 
minimal admixture of pure noise to the maximally entangled state such that a local realistic 
description is still possible. Operational realisation of the two quiVit measurement exhibiting 
strong violations of local realism involves entangled photons and unbiased multiport beam- 
splitters. The approach, extending at present to = 9, neither involves any simplifications, 
or additional assumptions, nor does it utilise any symmetries of the problem. 
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Chapter 1 
Introduction 



Bell's theorem |Pell64|| , formulated in 1964, initiated and revitalised several branches of 
modern physics. The paper was the first one to show that quantum entanglement cannot 
be in any way simulated by classical correlations. Within few years, a new branch of exper- 
imental physics emerged: multi-particle quantum interferometry. Since then it has evolved 
and extended its field of interest from two-photon to multi-photon correlations. Recently 
Bell-EPR correlations were observed for entangled atoms |[Hagley97| . For as much as nearly 
20-25 years the paper of Bell was studied mainly by people interested in the foundational- 
interpretational problems of quantum theory. Suddenly with the discovery of the possibility 
of employing Bell-EPR correlations | Eckert91 | in "quantum cryptography" ||Bennett84 ] and 
with the realisation of the importance of entanglement in the hypothetical quantum com- 
puters [[Feynman82| , peutsch89| , it turned out that the paper of Bell can be thought of as 
the first one in the field of quantum information. 

Studies of quantum information led to a proposal, employing entanglement, of quan- 
tum teleportation [ [Bennett^ . This phenomenon was observed in 1997 and seems to be 
at the moment the crowning achievement of quantum interferometry [Pouwmeester97|| . In- 
terestingly, the method to obtain quantum teleportation of photon's polarisation was de- 
veloped independently, as a by product of theoretical and experimental research towards 
obtaining Bell-EPR phenomena for particles originating from independent sources [Yurke91^ , 
Zukowski93a| , [Pan98|| . 

The same method was applied to obtain the first ever observations of Greenberger-Horne- 
Zeilinger correlations (GHZ) ||Bouwmeester99 1 . The 1989 theoretical discovery of GHZ cor- 



relations ||Greenberger89|| and the drastic amplification of the Bell theorem, which is implied 
by them, was the event in the research into foundations of quantum theory which am- 
plified interest in entanglement (with the interesting sociological consequence: the earlier 
terminology- correlated state- was replaced by the Schrodinger's term entanglement). 

The technological progress of 1980's and 1990's has enabled all that experimental and the- 
oretical activity to fiourish. The phenomenon of parametric down conversion (PDC) turned 
out to be a versatile source of entangled photons. The simplicity of the phenomenon of PDC 
has lead to an explosion of the number of experiments studying various aspects of entangle- 
ment or the basic phenomena linked with quantum information and quantum communication 
(dense coding | Mattle96 |, Bell-state measurement | Michler9^] , etc.). Pulsed down conver- 
sion enabled to observe two-photon interferometric effect for independently emitted photons 
Pan98|| . The future of experimental quantum information most probably will be associated 



with trapped atoms and microcavity-atom interactions ||Hagley97[| - fields in which extensive 
studies of entanglement are currently carried out. 

Interestingly, all these developments led to studies concerning entanglement properties 
of mixed states. In the case of the Bell theorem the pioneering works were by Werner 
(1989) |[Werner8"9| and the Horodecki Family (1995) ||Horodecki95|] . In recent years one 
observes an avalanche of works on the problem of separability of density matrices initialised 
by Peres |[Peres96|| and again the Horodecki Family (see e.g. ||Horodecki96|] ). The research 
into the separability has shown once more (recall Bell-Kochen-Specker theorem | Kochen6^] ) 
the quahtative difference between systems described by 2-dimensional Hilbert space (qubits) 
and those described by Hilbert spaces of higher dimension (quMts). 

The present work tries to answer some questions on the relation of the Bell theorem with 
various performed or proposed multiparticle (essentially, multiphoton) quantum interference 
experiments. 

The first part deals with some of performed experiments. Proposals of improvements and 
re-interpretations are given |[l|, |^ . 

Next, in part two, we study new methods of deriving Bell inequalities both for the 
standard two qubit experiments as well as to multi-qubit GHZ experiments 0, |]. Derivation 
of GHZ paradoxes for gedanken experiments involving M entangled quA^its observed beyond 
multiport beamsplitters ||Zukowski97b| j is also shown p, |[. 

Finally the last part is devoted to the Bell theorem without inequalities via a numerical 
approach utilising linear optimisation R, 0] . 
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Chapter 2 

Some history and basic notions 



2.1 Preliminaries 

According to a prevailing common opinion quantum mechanics is a fundamental theory which 
applies to all physical systems. Its predictive power is astonishing. Up to the present day 
there has not been a single experiment which invalidates it. However, the conclusions that 
can be drawn from quantum mechanics force us to entirely abandon the picture of nature 
implied by classical physics and the common sense. One of the main sources of difference 
between the quantum world and the classical one, in which we are particularly interested in 
this work, is entanglement. 

The notion of entanglement was introduced for the first time by Schrodinger to describe 
a situation in which 

Maximal knowledge of a total system does not necessarily include total knowledge 
of all its parts, not even when these are fully separated from each other and at 
the moment are not influencing each other at all... 



The work of Schroedinger [pchroedinger35|| was partially motivated by the seminal paper of 



Einstein, Podolsky and Rosen (EPR) ||Einstein"35| in which the authors used an entangled 



state (an EPR state) of two qubits to show that quantum mechanics could not be considered 
as a complete physical theory aiming to describe the phenomena occurring in micro world. 
Although they did not state it clearly, EPR effectively postulated the existence of local hidden 
variables, in the form of "elements of reality", which were to play the same role in quantum 
mechanics as the positions and velocities of particles in statistical classical mechanics and 
that were to solve the interpretational problems of quantum mechanics. That paper directly 
influenced the formulation of the Bell theorem^ Pell64| , Pell66| , |Bell87| 



In his 1964 paper Bell for the very first time showed^ that the idea of local hidden 
variables was in contradiction with quantum mechanics and, what is even more important, 
that it could be tested experimentally! That way the subject mainly discussed by physicists 
at the parties was brought to the realm of experimentally verifiable physics. 



"'^ According to Stapp [3tapp77|, the Bell theorem is one of the most important discoveries in modern 
physics. 



theorem by von Neumann [ vonNeumann32 1 excluding the possibility of the existence of hidden variables 
was formulated in 1930's but, as pointed out by Bell, the assumptions used by von Neumann were much too 
restrictive. 



2.2 Entanglement 



To discuss basic features of entanglement let us consider the following state of two two-state 
systems (qubits)^ 

IV') = ^(|0)i® 11)2 10)2), (2.1) 

where ® denotes the tensor product^ and kets |0)i, |l)i describe two orthogonal states of the 
i-th qubit. The above pure state describes a coherent superposition of two product states 
that occur with equal probability. 

According to quantum mechanics contains all available information about the state 
of the qubits. If we write ( |2.1| ) in the form of a density matrix pi2 

P12 = = \ [|o)ii(o| ® 11)22(11 - |o)ii(i| ® 11)22(01 

-|l)ii(0| ® 10)22(11 + |l)ii(l| ® 10)22(01] , (2.2) 
and trace out one qubit we obtain a density matrix describing the other qubit, which reads 

Pfc = ^(|l)fcfc(l| + |0)fcfc(0|), (2.3) 

A; = 1, 2. Such a density matrix describes a situation in which we have a chaotic mixture of 
two orthogonal pure states. This is a situation typical for the entanglement. We have the 
full possible information about the state of two qubits as a whole but we do not have any 
information about the state of each qubit separately, which in fact is not even defined! In 
addition, the properties of the qubits are tightly correlated. 

To see this more clearly let us consider the measurement of two dichotomic observables 
with spectrum consisting of ±1, which spectral decomposition has the following form 



with /c = 1, 2 and 



Ok{,<Pk) = |+, 4>k)kk{+, 4>k\ - |-, 4>k)kk{-, 4>k\ (2.4) 



\+Ak)k = ^{\0)k + e''^''\l) 



\-Ak)k = ^{\0)k-e'^''\l)k), (2.5) 

where 0^ G [0,27r]. The mean value Eqm{4>i,4'2) of the joint measurement of the observ- 
able Oi(0i) on the first qubit and the observable 02(^2) on the second one- the so called 
correlation function- reads 

^Qm(01,02) = (V^|Ol(0l)O2(02)|V^) = -COS(01 +02). (2.6) 



^If we take as the qubits two spin ^ particles and put |0) = | — |1) = | + ^) the corresponding state is 
a rotationally invariant state with the total spin equal zero- the so called singlet state. 
^This symbol will be only used when it makes the notation easier to read. 
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It is easy to see that, for example, whenever the sum of the phases (pi and 02 is modulo 
27r we observe the so called perfect correlations between the results of the measurements 
performed by the observers; if the first observer obtains +1 as the result of his measurement 
the second one obtains —1 and vice versa. However, each observer alone measures +1 and 



-1 with equal probability, which can be easily seen using the density matrix (|2.3| ) 



{Ok{<pk))p, = Tr{Ok{<pk)pk) = 0. (2.7) 



The above formulas clearly demonstrate that all information about the state is 
contained in the joint properties of the qubits. 



2.3 Elements of reality 



In their famous paper | Einstein"35[| EPR used the perfect correlations observed when measur- 
ing local observables on an entangled system of spatially separated qubits (|2.1|) to demon- 
strate that quantum mechanics is incomplete, i.e., that one needs some additional parameters 
to fully describe phenomena occurring in micro world. We briefly present their reasoning in 
the version of Bohm ||Bohm5l] , Pohm52|| to introduce the notion of local realism, which plays 
a central part in the Bell theorem. 

EPR reasoning goes as follows. According to quantum mechanics all we know about the 
system of two entangled qubits is encoded in the state (|2.1|). Quantum mechanics also tells us 
that we cannot consider simultaneous measurements of two non commuting observables and 
therefore does not even define predictions for such cases. Let us imagine that the observers 
one and two are spatially separated and that they simultaneously measure the observables 
Oi(0i = 0) and 02(^2 = 0) (see the picture (^]l|)) on subsystems described by the singlet 
state. 



A2 
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Figure 2.1: The scheme of the Bell type experiment. Solid arrows represent actual positions 
of the measuring apparatus, dashed arrow the future ones. Detectors Ai,A2,Bi, B2 register 
incoming qubits. 



For such observables the perfect correlations occur, which means that if the first observer 
has obtained +1 the second one, because of ( p.6|) , must have obtained —1 and vice versa. 
Moreover, due to the spatial separation of the observers and the non superluminal velocity 
of propagation of any interaction (information) in nature (locality), the outcome of the 
first (second) observer cannot be influenced by the choice of observable measured by the 
second (first) observer and neither by its outcome. All correlations between the results of 
measurements must have been established in the source. 
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However, the second observer could have measured the observable 02(02 — f ) instead of 
the previous one characterised by 02 = 0. Then, the result of this would-be measurement, 
would have enabled him to infer with probability equal to one and without disturbing the first 
qubit the result of the measurement of the observable Oi(0'^ = ^) by the first observer. At 
this stage of reasoning EPR introduce the notion of physical reality pinstein35| : 



If, without in any way disturbing a system, we can predict with certainty (i.e. 
with probability equal to unity) the value of a physical quantity, then there exists 
an element of physical reality corresponding to this physical quantity. 

According to this definition the inference made by the second observer^ about the result of 
the measurement that could have been made by the first observer has a well defined physical 
meaning. This would mean that it is possible to ascribe the definite values to the results of 
the measurement of two non commuting observables Oi(0j = 0) and Oi(0'j = |) (i = 1,2), 
the statement which makes quantum purist's hair stand on ends. Therefore, according to 
EPR, quantum mechanics is incomplete and should be completed at least by introducing 
elements of reality into the description of quantum state. 



2.4 The Bell theorem 

The hypothesis of local hidden variables, which effectively stems from the notion of elements 
of reality, was proved to be unacceptable for systems described by quantum mechanics by 



Bell, [[Bell64|] who found an inequality which should be obeyed by any local hidden variable 
theory, and which is violated by predictions of quantum mechanics. Here we present the 
derivation of the Clauser-Horne-Shimony-Holt (CHSH) inequality ||Clauser69|] . 



To this end, let us again consider the experiment in which two spatially separated ob- 
servers perform the measurements of the observables defined in ( p.4|) on the state (|2.1| ). Each 
of them as a result of their measurement obtains only one of the two possible values ±1 (they 
measure bivalued observables). Although the different results of the measurement at each 
observer's side occur with equal probability the results of joint measurements are correlated, 
which is expressed by the formula ( p.6|) . If deterministic local hidden variables exist they 
predetermine the results of each single measurement for each observer at the time of the 
emission of each pair of qubits. In other words, the value of the local hidden variable for a 
particular emission of a pair of qubits would allow us to predict the result of the measure- 
ment made by each observer with certainty. However this value is hidden. To model the 
probabilistic nature of quantum experiments we assume that there exists some probability 
distribution of local hidden variables associated with a given quantum mechanical state, 
which represents our lack of knowledge about them. 

To express the above idea of local hidden variables mathematically we assume that there 
is a set of hidden variables A on which we can define a probabilistic measure dp{X). We also 
assume the existence of two bivalued functions Oi, O2 defined on the space A which take only 
the values ±1. Each of these two functions must depend on the local parameters 0i and 02 
characterising the experiment performed by the first and the second observer respectively. 
Because of the assumption of locality the function Oi can solely depend on the 0i and the 



^The reasoning can be reversed and the inference can be made by the first observer about the second one. 
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other one on 02- Thus, the correlation function based on the idea of local hidden variables 
must have the following form 



^w(0i,02)= / rfp(A)Oi(A,0i)O2(A,02), (2.8) 

J A 

where Jj^dp{X) = 1. 

Now, let us imagine that each observer performs two mutually exclusive experiments 
characterised by 0j, 0'j and let us consider the following expression made out of the four local 
hidden variables correlation functions 



Clhv = Elhv{(I)i, 02) + Elhv{(P'i, 02) + ELHvih, 02) - Elhv{<I)'i, 02) 

dp{\) [Oi (A, 0i)(O2(A, 02) + 02(A, 0'2)) + Oi(A, 0;) (02(A, 02) - 02(A, 0'2))](2.9) 



It is easy to see that the modulus of the expression in the square brackets is either —2 or 
+2. Therefore, the hypothesis of local hidden variables implies that the following inequality 
(Bell-CHSH inequality) must be valid 

-2 < Clhv < 2. (2.10) 



Is the CHSH inequality always satisfied by quantum predictions for (pTT|)? To answer 
this question let us put 0i = 0, (p[ = | and 02 = — 7r/4, 02 = 7r/4. For these values of local 
parameters one has 

Cqm = Eqm{4>1, 02) + ^QA/(0'l, 02) + ^QA/(01, 02) " ^QAf (0'l, 02) 

= - COS(01 + 02) - COS(0; + 02) - cos(0i + (f)^) + cos{(f)[ + (p'^) = -2^2. (2.11) 

Because —2^/2 < —2, we have a contradiction. 

The above result, known as the Bell theorem^, needs some further explanation. In our 
reasoning we have made two crucial assumptions without which the theorem would not be 
valid. These assumptions are: locality and realism. The Bell theorem tells us that either 
notion of locality, or realism, or both are false in quantum theory ||Redhead87| . 



Another remark is that the Bell-CHSH inequality can be directly applied to any ex- 
perimental data. Also, even if quantum mechanics is not valid and we will find another 
better theory we can still, using the Bell inequality, verify whether this new theory fulfils 
the necessary condition for the local realistic description or not. 

The final remark is that one may consider the existence of the so called stochastic local 
hidden variables [|Clauser74|| , which do not predict with certainty the results of local mea- 



surements but give merely the probabilities of their occurrence. In such a case instead of 
functions 0„ (n = 1, 2) appearing in ( |2.8| ) we have the probabilities Pi(m|A, 0i), P2("^'|A, 02) 
giving the ratio of occurrence of the results m, m' when measuring observables characterised 
by parameters 0i, 02 respectively (obviously they sum up to one, i.e., J2m=±i -Pi('^|A, 0i) = 1, 
i = 1,2). The relation between Pj(m|A,0i) and Oi(A,0j) {i = 1,2) is such that within this 
description Oj = ±1 have to be replaced by 

O,(A,0,)= J2 mPiim\X,(l)i). (2.12) 

m=±l 



®The inequality which must be obeyed by any local and realistic theory is usually called the Bell inequality 
whereas the violation of such an inequality by quantum mechanics is called the Bell theorem. 
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with the values of modulus bounded by 1. With such probabilities the idea of stochastic 
local hidden variables is to reproduce the quantum probabilities PQM{''Ti,m'\(f)i,(f)2), i-e., the 
probabilities of obtaining the result m and m' by the first and the second observer when 
measuring the observable characterised by 0i and 02 respectively, by local hidden variables 
probabilities of the form 



PHv{m,m'\4>i,(t>2) = J dp{X)Pl{m\X,(t>l)P2{m'\X,(f)2)■ 



{2.13) 



It is clear that any deterministic local hidden variables theory can be always treated as a 
stochastic one. Fine [[Fine82|| proved that a stochastic local hidden variable theory implies 
the existence of an underlying deterministic oneQ 

In general, all Bell-type inequalities found since the famous Bell paper ||Bell64|| constitute 
only necessary conditions for the existence of local hidden variables. The exceptions are the 
full set of four Clauser-Horne inequalities (CH) ||Clauser7^ and the Bell-CHSH inequality, 
which were proved by Fine |[Fine82|| to be also sufficient ones for dichotomic observable^ 
(see also [[PeresOQII ). 

In 1989 Greenberger, Horne and Zeilinger (GHZ) ||Greenberger89|| used a maximally en- 
tangled state of four qubits to show that the discrepancy between local realism and quantum 
mechanics is much stronger than that observed in two qubit correlations|. By a clever trick 
they showed that the idea of local realism breaks already at the stage of defining elements 
of reality. 

The Bell theorem does not have to be restricted to two or three qubit correlations. The 
discrepancy between local realism and quantum mechanics can be also proved for entangled 
particles each described by an dimensional Hilbert space- so called quA^its (see, for in- 
stance, [[Mermin8CI| , [MerminS^ , |Garg82| , pisin92| , [Peres92| , |W6dkiewicz94| ] )- as well as for M 
entangled qubits (see, for instance, [[Mermin90b|| ) . 



2.5 Experimental tests of the Bell theorem 
2.5.1 First experiments 

The first experimental test of Bell inequality was performed by Freedman and Clauser 
[b'reedman72|| with photons from atomic cascade decays. They observed violation of Bell 
inequality and confirmation of quantum mechanical predictions. However, in the experiment 
the detection efficiency and the angular correlation of the photon pairs was low, and no care 
was taken to make the two polarising settings detection stations to be set independently 



Clauser 74, 5antos92 



Such care was taken in the most quoted experiment by Aspect, Grangier and Dalibard 
Aspect8^] . In this experiment fast switchings of the analyser position to prevent "commu- 
nication" between the source and the analyser was used. However, as it was pointed out in 
[Zeilinger86|| the periodic switching was not truly random and was predictable after a few 



periods of the switchQ. 



''We do not take into account non Kolmogorovian probability calculus. 
®The CHSH inequality is sufficient if one assumes certain symmetries of the probabilities. 
^Later Mermin simplified the proof and derived GHZ paradox for three qubits Mermin90a| . 
^"in 1998 Weihs et. al. \ Weihs9^ for the first time performed an experiment in which true locality 
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In all experiments (except one, with systematic errors) the violation of the Bell inequal- 
ities (with certain additional assumptions) was observed. An excellent review concerning 
these pioneering experiments can be found in [ClauserTS . 



2.5.2 New experiments 



Recently the experiments, in which the entangled pairs of photons are generated in the 
process of parametric down conversion (PDC), dominated the field of laboratory tests of 
local realism. In the PDC process the pairs of photons are spontaneously created. The 
propagation directions and the frequencies of created photons (the photons are called idlers 
and signals) are highly correlated, which is used to generate an entangled state. In the 
type-II down conversion one has also correlated polarisations. 

Among the Bell-type experiments with PDC process the following ones will be mainly 
addressed to in this work 



Alley-Shih and Ou-Mandel experiments 



UvM 



entanglement swapping | [Pan98 | 



GHZ experiment [ |Bouwmeester99| | . 



First Bell-type experiments with a PDC source of correlated photons were experiments 
by Alley, Shih and Ou, Mandel. In both of them the violation of the Bell inequality (by 
3 standard deviations in ||Shih88| ] and 6 standard deviations in ||Ou88|| ) and confirmation of 
quantum mechanics was reported. However, in both experiments only coincident counts were 
measured (half of the events were discarded). That raised some doubts about the validity of 
the experiments as tests of local realism [^antos92 , Caruccio94 |. The situation was clarified 
in |Popescu^7|] where it was shown that one does not need to discard "wrong" events to test 
local realism in the experiments. In this dissertation we perform further theoretical analysis 
of these experiments (see Part I). 

In 1993 Zukowski et. al ||Zukowski93a|| showed experimental conditions to entangle par- 
ticles (photons) originating from independent sources]^ Five years later, in 1998, the first 
entanglement swapping experiment was performed ||Pan98| | . The visibility of around 65% 
was observed (the notion of visibility is explained in the next subsection). 

In 1999 Bouwmeester et. al | [Bouwmeester99| ] reported the first experimental observation 
of the GHZ correlations. The experiment was based on the techniques developed in the tele- 



portation experiment | Bouwmeester97 | and entanglement swapping experiment |Pan98 ]. In 
the experiment pairs of entangled photons (entangled polarisations) produced in a nonlinear 
crystal pumped by a short pulse of ultraviolet light from the laser were used. The applied 
technique to obtain GHZ correlations rests upon an observation that when a single particle 
from two independent entangled pairs is detected in a manner such that it is impossible to 
determine from which pair the single came, the remaining three particles become entangled. 
The high visibility of around 60% was observed. 



condition was enforced. In the experiment two observers were spatially separated by the distance of 400in, 
which means that the time of the measurement had to be shorter than 1.3/Lts to prevent communication with 
the speed of light between observers. They succeeded to achieve the time of measurement within 100 ns and 
the violation of the CHSH inequality by 30 standard deviations was observed. 
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The first proposal was given by Yurke and Stoler |Yurke92 
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The experimental realisation of entanglement does not have to be restricted to massless 
particles (photons). Hagley et. al. reported the experiment in which entangled atoms were 
produced ||Hagley97[] . They demonstrated the entanglement of pairs of atoms at centimetric 



distances and measured their correlations. The visibility of only around 25% was observed. 

2.5.3 Problems encountered in the Bell type experiments 

However, in all experiments performed thus far there has been the problem with a low quan- 
tum efficiency rj of detectors used to register incoming particles^. The quantum efficiency is 
defined as the ratio of the number of detected particles to the number of the emitted ones0. 
If it lies below the threshold value Vj^ = 2a/2 — 2 there is no violation of the Bell-CHSH and 
the CH inequality for maximally entangled two qubits. 

Since the collection efficiency in all experiments done so far was much lower than rjl^, 
all claims about the violations of Bell inequalities in performed experiments are based on 
the assumption that the observed sub ensemble of particles is a representative one for the 
emitted ensemble (so called "fair sampling assumption"). 

In the real experiment one usually cannot obtain a pure maximally entangled state of 
two qubits due to some imperfections in the source producing the state and other difficulties. 
As a simple generic model of such experimental imperfections one can take 

p{F2) = F2Pnorse + (1 " i^2)Pp«re, (2.14) 

where for qubits Pnoise = {I is a unit 2x2 matrix), Ppure = (fo^' the definition of 

see ( p.l| )) and the real parameter F2 (the index 2 stands for qubit) lies between zero and 
one (0 < F2 < 1). The pnoise describes a totally chaotic mixture of two qubits. Results of any 
measurements carried out on the pnoise are completely uncorrelated therefore the parameter 
F2 can be interpreted as the number telling us how much noise is contained in the system. 
It is easy to check that for p{F2) the correlation function defined in ( p.6| ) reads 



^QAf (01, 02) = Tr[p(F2)6i(0i)62(02)] = -(1 - F2) cos(0i + 02). (2.15) 

We see that if F2 > the amphtude of the correlation function is reduced. 

In quantum interferometry the number 1 — F2 is directly linked with the visibility (con- 
trast) of interferometric two-qubit fringes]^. 

Therefore, sometimes it is convenient to consider the parameter V2 = 1 — -^2 instead of 
F2, which is usually done in the description of real experiments. Throughout the dissertation 
both parameters will be used depending on the context. 

^^More precisely the quantum efficiency describes the full detection stations, including all devices that 
collect the incoming radiation (lenses, etc.). 

-'^^In the operational terms it is defined for a detection station A as the ratio of coincident correlated counts 
at the pair of detection stations A and B to the ratio of singles at the station B. 

^''The traditional definition of interference visibility (contrast) is given by 

V = - —j— l^-iOJ 

^max r -'-min 

where / is the intensity (of light) in an interference pattern. Imax refers in the case of spatial pattern to the 
maximum intensity and Imin to its neighbouring minimum. The same formula can be used for an output of a 
Mach-Zehnder interferometer into one of its exit arms. In this case Imax is the maximum intensity and Lm-in 
its neighbouring minimum which occurs after suitable change of the phase shift. In the case of a quantum 
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In the entanglement swapping experiment |Pan98|| and in the experiments with atoms 
H.agiey97| ] the observed visibihty was quite low: around 65% in the entanglement swapping 
and 25% in the experiment with entangled atoms. 

It is obvious that if there is too much noise in the system, i.e., the visibility is low, then 
one cannot observe violations of Bell inequalities. 

To summarise, if in a Bell-type experiment with two maximally entangled qubits V2 < ^ 

{F2 > ^~2^) the Bell-CHSH inequality (and also the CH inequality) cannot be violated. 
Furthermore, in the case of the Bell-CHSH inequality, if r]2 < 2^ - 2, then even if we have 
the perfect case, i.e., V2 = 1 {F2 = 0), the violation of the inequality has only a bona fide 
status- one has to use the fair sampling assumption. In all experiments thus far performed 
one has r]2 <^ 2\f2 — 2. Similar situation (to be described in more details later) occurs for 
the GHZ correlations (with new specific threshold r^'s and V^'s). 



process for a single particle the visibility is defined, in analogy, as 



Prr 



P 

J rr. 



Prr 



Prr. 



(2.17) 



where Pmax and Pmin ^I'e the maximal and minimal probabilities for detection of the particle in a specified 
output of an interferometer (also obtainable for certain different phase settings). 



The visibility of two particle interference is again defined using the same general rule (2.17). However, in 
this case P's refer to the probability of coincident counts at a pair of detectors. In the case when both single 
particle and two particle interference occurs in the experiment, the relation between the two visibilities is 
quite subtle, and the two-particle visibility has to be redefined |Jaeger93 . However, throughout the present 



work we shall discuss only the cases for which no single particle interference occurs. 
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Part I 

New theoretical analysis of 
AUey-Shih, Ou-Mandel and 
entanglement swapping experiments 
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* 



Some of the performed Bell-type experiments have a distinguishing trait. Not all event 
observed follow the standard pattern assumed in the usual derivations of Bell inequalities. 
To such experiments belong the AUey-Shih-Ou-Mandel experiment ||Shih88| , |Ou88| 
entanglement swapping experiment [[Pan98| . 



and the 



In all these experiments even in the ideal, gedanken case, half or more of the emissions 
do not lead to correlated counts at spatially separated detector stations. Therefore in order 
to prove that such experiments can be indeed considered as tests of local realism one has 
to perform an analysis which takes into account this characteristic trait. Such an analysis 
will be presented below for the Alley-Shih, Ou-Mandel experiments and the entanglement 
swapping experiment]^. Special care will be taken to discuss the question of whether two 
versus single photon counts distinguishability is the necessary requirement for the studied 
experiments. 



"'^^The experiments proposed by Franson (1989) share properties which were thought to be of a similar 
nature. By a closer inspection it turns out that they are different. For explanation see | Aerts99[ . 
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Chapter 3 

AUey-Shih and Ou-Mandel 
experiments: resolution of the 
problem of distinguishability of single 
and two photon events [1] 



3.1 Introduction 

The first Bell-type experiments which employed parametric down conversion process as the 
source of entangled photons were those reported in refs ||Shih8^] and ||Ou88|| . However, the 



specific traits of those experiments have led to a protracted dispute on their validity as tests 
of local realism. In this case the issue was not the standard problem of detection efficiency 
(which up till now permits a local realistic interpretation of all performed experiments). 
The trait that distinguishes the experiments is that even in the perfect gedanken situation 
(which assumes perfect detection) only in 50% of the detection events each observer receives 
a photon, in the other 50% of events one observer receives both photons of a pair while the 
other observer receives none. The early "pragmatic" approach was to discuss only the events 
of the first type (as only such ones lead to spatially separated coincidences). Only those were 
used as the data input to the Bell inequalities in |^hih88|| and |Pu88|| . This procedure was 



soon challenged (see e.g. [|Kwiat94| , |Kwiat95[] , and especially in the theoretical analysis of ref. 



CjaruccioQ^ ) , as it raises justified doubts whether such experiments could be ever genuine 



tests of local realism (as the effective overall collection efficiency of the photon pairs, 50% 
in the gedanken case, is much below what is usually required for tests of local realism). Ten 
years after the first experiments of this type were made, finally the dispute was resolved 
Popescu97[] . It was proposed, to take into account also those "unfavourable" cases and to 



analyse the entire pattern of events. In this way one can indeed show that the experiments 
are true test of local realism (namely, that the CHSH inequalities are violated by quantum 
predictions for the idealised case). The idea was based upon a specific value assignment for 
the "wrong events". However, the scheme presented by Popescu et al ||Popescu^7| has one 



drawback. The authors assumed in their analysis that the detecting scheme employed in 
the experiment should be able to distinguish between two and one photon detections. This 
was not the case in the actual experiments. The aim of this chapter is to show that even 



on 



this is unnecessary, all one needs is the use of the specific value assignment procedure of 
[FopesculJ7[ . 

Finally, we shall also give prediction of all effects occurring in the experiment. It is quite 
often overlooked that a kind of Hong-Ou-Mandel dip phenomenon ||Hong87| | can be observed 
in the experiment. 



3.2 Description of the experiment 

In the class of experiments we consider (see [Popescu97|| a type I parametric down- 



conversion source 



Hong85[| is used to generate pairs of photons which are degenerated in 



frequency and polarisation (say x) but propagate in two different directions. One of the 
photons passes through a wave plate (WP) which rotates its polarisation by 90". The two 
photons are then directed onto the two input ports of a (nonpolarising) "50—50" beamsplitter 
(BS). The observation stations are located in the exit beams of the beamsplitter. Each 
local observer is equipped with a polarising beamsplitter^, orientated along an arbitrary axis 
(which, in principle can be randomly chosen, in the delayed-choice manner, just before the 
photons are supposed to arrive). Behind each polarising beamsplitter are two detectors, Df, 
Di and D2, D2 respectively, where the lower index indicates the corresponding observer 
and the upper index the two exit ports of the polarised beamsplitter ( " + " meaning parallel 
with the polarisation axis of the beamsplitter and " — " meaning orthogonal to this axis). 
All optical paths are assumed to be equal. 



3.3 Quantum predictions 

Let us calculate the quantum predictions for the experiment. We will use the second quan- 
tisation formalism, which is very convenient here, since the whole phenomenon observed in 
the experiment rests upon the indistinguishability of photons. 

After the action of the wave-plate one can approximate quantum mechanical state de- 
scribing two photons emerging from a non - linear crystal along the "signal" and the "idler" 
beam by 

l^o) = «UL^|0), (3.1) 

where a\g and a\g are creation operators and |0) denotes the vacuum state. Subscripts x,y 
decode the polarisation of the photon (either along x or y axis). The beamsplitter action 
can be described by 



73(4 + 4)' 



(3.2) 



^Following references |Garuccio94| and |Popescu97| we assume that both local detection stations are 
equipped with polarising beamsplitters, and each of the output ports is observed by a detector. In the actual 
experiments [3hih88| and |Ou88| at each station only one of the outputs was monitored. 
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Figure 3.1: Schematic of the setup. For explanations see the main text. 

where ct, dt, ct, rft are operators describing output modes of the beamsphtter (c stands for 
the first observer and d for the second one). Thus our state |^o) changes to : 

1^) = ^(44-44+44+^44)10) (3-3) 

Next comes the action of the polarisers in both beams, which can be described as 

nt = cos(^i)n|| + sin(^i)nj^ 
nt = sin(6'i)n|| - cos{9i)'n}j_, 

(3.4) 

where = or ci^, and describes the mode parallel to polarizer's axis and describes 
the mode perpendicular to polarizer's axis; 6 is the angle between the x axis and polarizer's 
axis. Thus the final state reaching the detector reads 

iV'/mai) = |[sin(6'i - e2)\c\\,d\\) 
+ cos(6'i - 92)\c\\di_) 
- cos(^i - ^2)|c±, d\\) + sin(^i - ^2)|c±, c?_l) 
+2^sin(2^i)|2c||) +i^sin(2^i)|2ci) 
-i cos(2^i)|cx, C||) + i-^ sin(2^2)|2(i||) 



oo 



+i^sin(2^2)|2rfx) -«cos(2^2)M||,rf±)], 

(3.5) 

where e.g. |c||, d\\) denotes one photon in the mode cy, and one in rfy, whereas |2c||) = ^C|| |0) 
denotes two photons in the mode cy . 

Let us denote by P{i, 9i;j, 62) the joint probabihty for the outcome i to be registered by 
observer 1 when her polariser is oriented along the direction that makes an angle 9i with the 
X direction and the outcome j to be registered by observer 2 when her polariser is oriented 
along the direction that makes an angle 62 with the x direction. Here z, j = 1 — 6 and have 
the following meaning [[Popescu^j : 



l=one photon in no photon m. 
2=one photon in , no photon in D~ 
3=no photons 

4=one photon in and one photon in D~ 
5=two photons in no photon in D~ 
6=two photons in D~ , no photons in . 

The quantum predictions for joint probabilities of those events are given by: 



P(l,^i;l,^2) = P(2,0i;2,02) 
P(2,^^i;l,^2) = P(l,^i;2,e2) 



cos 2(01-^^2)], 

|[l + cos2(^^i-^^2)], 



P(5, 3, 62) = P(6, Or, 3, 62) = \ ^^^26^), 
P{3, 5, 62) = P(3, ^i; 6, 62) = \ sm\2d2), 
P(4,^^i;3,^^2) = icos2(2ei), 
P(3,^i;4,^2) = icos2(2e2). 



(3.6) 
(3.7) 
(3.8) 
(3.9) 
(3.10) 
(3.11) 



Following ||Popescu^7| we associate with each outcome registered by the observer 1 and 2 a 
corresponding value Oj and hj respectively, where ai = 61 = — 1 while all the other values are 
equal to 1. Let us denote by E{Oi, O2) the expectation value of their product 



E{9^,e2) = Y.a,hjP{i,er-j,e2) 

After simple calculations one has: 



(3.12) 



^(V'l,V'2) 

= -icos(V^i + V^2) + i (3.13) 

where we have put 29k = (— l)'^~^V'A:• 
The above formula for the correlation function is valid if one assumes that it is possible to 
distinguish between single and double photon detection. This is usually not the case. Thus 
it is convenient to have a parameter a that measures the distinguishability of the double and 
single detection at one detector ( < a < 1, and gives the probability of distinguishing by 
the employed detecting scheme of the double counts). The partial distinguishability blurs 
the distinction between events 1 and 6 (2 and 5) and thus part of the events of the type 6 are 



interpreted as of type 1 and are ascribed by the local observer a wrong value, e.g. an event 
of type 6, if both photons go to the " — " exit of the polariser, can be interpreted as a firing 
due to a single photon and is ascribed a —1 value. Please note that such events like 1 or 2 
in station 1 accompanied by 3 (no photon) at station 2 do not contribute to the correlation 
function because for any a P(l, 6'i; 3, 62) = P{2, 9i; 3, 62). 

If the parameter a is taken into account the correlation function acquires the following 
form: 



E{ilJi,ilj2;a) 

■\ cos('?/'i + ll>2 

~2 , 



1, 



+ - a)(cos2 + cos^ ^2 



(3.14) 



3.4 Conditions to violate local realism 

After the insertion of the quantum correlation function (|3.14|) into the CHSH inequality, 



-2 < E(7/;i,V'2;a) + ^(V'l,^2;a) 
+E(^i,V'^;a)-E(^l,V^^;a)) <2, 

one obtains: 

-2 < -i[cos(V^i + ^2) + cos(V^; + ^2) 

+ COs('?/'i + ■j/'a) - COs('?/'i + ^'2)\ + OL 

+ |(l-a)(cos2^/>i + cosV2) <2. (3.15) 

The interesting feature of this inequality is that it can be violated for all values of a. 
What is perhaps even more important, it can be robustly violated even when one is not 
able to distinguish between single and double clicks at all (a = 0). The actual value of the 
CHSH expression can reach in this case 2.33712 (a numerical result), which is only slightly 
less than the maximal value for a = 1, which is a/2 + 1 fa 2.41421. Therefore we conclude 
that in the experiment one can observe violations of local realism even if one is not able 
to distinguish between the double and single counts at one detector. That is, the essential 
feature of the method of [[Popescu^j to reveal violations of local realism in the experiment 
of this type is the specific value assignment scheme and not the double-single photon counts 
dist inguishability. 

The specific angles at which the maximum violation of the CHSH inequality is achieved 
for a = differ very much from those for a = 1 (for which the standard result is reproduced), 
and they read (in radians) ipi = 2.93798, ip[ = 4.25513, = -0.20241 and ip2 = 1.11708. 

Let us notice that with the setup of ( p.l|) one is able to observe effects of similar nature to 
the famous Hong-Ou-Mandel dip ||Hong87 . These are revealed by the probabilities pertain- 



ing to the wrong events (3.8-3.11). Simply for certain orientations of the polarisers, if the two 
photons emerge on one side of the experiment only, then they must exit the polarising beam- 
splitter via a single output port (this effect is due to the bosonic-type indistinguishability of 
photons, see |[Hong87||). 
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Finally let us discuss what is the critical efficiency of the detection of the experiments 
of this type. To this end, in our calculations we will use a very simple model of imperfect 
detections: we insert a beamsplitter with reflectivity a/1 — t], in front of an ideal detector, 
which observes only the transmitted light. This results in the system behaving just like a 
detector of efficiency r]. If we assume that the incoming light is described by a creation 
operator then transmitted mode is denoted as t]j whereas reflected mode is denoted as 
and one has 

at = yT^rlt + v^tlt. (3.16) 
For instance, if one takes the following part of the state vector ( |3.5| ): 

cf|4||0). (3.17) 

the beamsplitter model of an imperfect detector transforms this term into: 

[(1 - vyi^d^^ + \jv{^-vy\^t\ 
+\/^(i-^)4„^dii+<,4|,]|o). 

(3.18) 

The probabilities now read: 



P(3,ei; 


2,^^2) 


= P(2,^^i;3,02) 






P(i,^i; 


3,^^2) 


= P(3,0i;l,02) = 


= r/(l -r]) 


(3.19) 


ni,^i; 


1,^2) 


= P(2,0i;2,^^2) = 




(3.20) 


P(2,ei; 


1,^^2) 


= P(1,0;2,02) = 




(3.21) 


P(5,ei; 


3,^^2) 


= P(6,0i;3,02) = 


- \v'[Mmf 


(3.22) 


P(3,^^i; 


5,02) 


= P(3,0i;6,^^2) = 


- ^r^^[sin(2^^2)]^ 


(3.23) 




3,^2) 


= \^'[cos{2e,)Y 




(3.24) 


P(3,^i; 


4,^2) 


= ir^^[cos(2^^2)]^ 




(3.25) 



The correlation function, which includes the inefficiency of the detection reads 

^(V'l, ^2; Tl, «) = V^E{ij,,ij2; «) + (1 - vY, (3.26) 

where £{1^1,1^2', a) is given by ( p.l4|) . We have tacitly assumed here that the parameters 
a and t] are independent of each other (this assumption may not hold for specific technical 
arrangements). Putting this prediction into CHSH inequality, assuming that a = 1 (full 
distinguishability) we obtain a minimum quantum efficiency needed for violation of local 
realism equal to 0.91, whereas for other values of a we have: for a = 77 = 0.926; for 
a = 0.5 rj = 0.92; for a = 0.75 rj = 0.92; for a = 0.875 77 = 0.91. One should note here 
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that the method of value assignment of [Popescu97|| is in accordance with the method given 
by Garg and Mermin | |Garg85 | for the optimal estimation of required detector quantum 
efficiency to violate local realism in a Bell-test. Thus the obtained efficiencies are indeed 
the lowest possible, and show that experiments of this type are not good candidates for a 
" loophole- free" Bell-test [^antosQ^ , nevertheless due to the fact that the whole observable 
effect is a consequence of quantum principle of particle indistinguishability such test are 
very interesting by themselves - they reveal the entanglement inherently associated with this 
principle. 



3.5 Conclusions 

To conclude, we state that the possibility of distinguishing between single and two photon 
detection events, usually not met in the actual experiments, is not a necessary requirement 
for the proof that the experiments of Shih- Alley and Ou-Mandel are, modulo fair sampling 
assumption, valid tests of local realism. We also show that some other interesting phenomena 
(involving bosonic type particle indistinguishability) can be observed during such tests. 



Chapter 4 

Better entanglement swapping [2] 



4.1 Introduction 



Until recent years it was commonly believed that particles producing EPR-Bell phenomena 
have to originate from a single source, or at least have to interact with each other. However, 
under very special conditions, by a suitable monitoring procedure of the emissions of the 
independent sources one can pre-select an ensemble of pairs of particles, which either reveal 
EPR-Bell correlations, or are in an entangled state. The first explicit proposal to use two 
independent sources of particles in a Bell test was given by Yurke and Stoler |[Yurke92 



However, they did not discuss the importance of very specific operational requirements nec- 
essary to implement such schemes in real experiments. Such conditions were studied in 
Zukowski93a|| and ||Zukowski95| . 



The method of entangling independently radiated photons, which share no common past, 
[ [Zukowski93a|| is essentially a pre-selection procedure. The selected registration acts of the 
idler photons define the ensemble which contains entangled signal photons (see next sec- 
tions). Surprisingly, such a procedure enables one to realize the Bell's idea of "event-ready" 
detection. This approach for many years was thought to be completely infeasible and thus 
no research was being done in that direction |piauser78| . This so-called entanglement swap- 
ping technique |[Zukowski93a|| , was also adopted to observe experimental quantum states 
teleportation ||Bouwmeester97 . 

High visi- 



The first entanglement swapping experiment was performed in 1998 [Pan98 



bility (around 65%) of two particle interference fringes were observed on a pre-selected subset 
of photons that never interacted. This is very close to the usual threshold visibility of two 
particle fringes to violate some Bell inequalities, which is 70.7%. Therefore there exists a 
strong temptation for breaking this limit, and in this way showing that the two particle 
fringes due to entanglement swapping have no local and realistic model. 

However, due to the spontaneous nature of the sources involved, the initial condition 
for entanglement swapping cannot be prepared. Simply the probability that the two sources 
would produce a pair of entangled states each is of the same order as the probability that one 
of them produces two entangled pairs. In the latter case no entanglement swapping results. 
Nevertheless, such events can excite the trigger detectors (which in the case of the right initial 
condition select the antisymmetric Bell state of the two independent idlers). Therefore they 
are an unavoidable feature of the experiment, and have to be taken into account in any 
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Figure 4.1: Entanglement swapping. Two type-II down conversion crystals are pumped by 
a pulsed laser. The radiation from each of the crystals is entangled in polarisations (e.g. if 
one has an H polarised photon in mode a, then in mode 6 is a V polarised photon). The idler 
photons (in modes d and e) are fed into a non-polarising beam splitter BS. Simultaneous 
firing of the trigger detectors Tg and pre-selects, in the event-ready way, a sub-ensemble 
of detection events behind the polarising beam splitters A and B. The orientation of the 
polarising beam splitters can be set at will by the local observers. The output signal photons 
are registered by the two local detection stations, consisting of detectors denoted by + and 



analysis of the possibility of finding a local realistic description for the experiment. 

The aim of this chapter is to perform such an analysis. We shall show that if all firings 
of the trigger detectors are accepted as pre-selecting the events for a Bell- type test^, one 
must necessarily, at least partially, be able to distinguish between two and single photon 
events at the detectors observing the signals to enable demonstrations of violations of local 
realism. Whereas, if one accepts additional selection at the trigger detectors, based on the 
polarisation of the idlers, detectors possessing this ability are unnecessary. We shall present 
our argumentation assuming that the reader knows the methods and results of [ Zukowski93a| , 



Zukowski95| , |Pan98|| . The analysis will be confined to the gedanken situation of perfect 



detection efficiency (the results can be easily generalised to the non-ideal case) 

4.2 Description of the experiment 



Consider the set-up of ( [4.1|) , which is in principle the scheme used in the Innsbruck experi- 
ment |[Fan98|| . Two pulsed type-II down conversion sources are emitting their radiation into 
the spatial propagation modes a and h (signals), c and d (idlers). Due to the statistical 
properties of the PDC radiation, the initial state that is fed to the interferometric set-up has 
the following form: 

i^) = E^=o(^(«l.4+44)r 



^As was the case in the actual experiment. 
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xE^=o(^(4&y + eUl/))"|0), (4.1) 

where, for instance, a}^ denotes the creation operator of the photon in beam a having "hori- 
zontal" polarisation. As for the entanglement swapping to work one cannot have too excessive 
pump powers [^ukowski99 |, the 7 coefficient can be assumed small. Therefore we select only 



those terms that are proportional to 7^, as these are the lowest order terms terms that can 
induce simultaneous firing of both trigger detectors. They read 

+ (al,4 + aU^H? + (eL&l- + eUL)')|0). (4.2) 

The factor ^7^ simply gives the order of magnitude of the probability of the two trigger 
detectors to fire, and therefore we drop it from further considerations. The action of the 
non-polarising beam sphtter (BS) is described by = -^{dl. + ie\) and e\ = :^(ei + idl) 

where x = H or x = V, and e and d represent the modes monitored by the trigger detectors 
behind the beam splitter. Taking into account only the terms in ( [4.2|) that lead to clicks at 
two trigger detectors we arrive at 

m = (^(ag + b^H)eUv + Kay + ^'y )44 

lia^^h^y - al.6t^)(et4 - et4))|0). (4.3) 
It is convenient to normalise and rewrite the above state into the form: 



+zV2{j=^al' + j=^bl')\HH) 

+ {^{aUl + b^nbl) - K^fel - 4b^ \HV), (4.4) 

where \VV) = 44 1 0), \HH) = 44 1 0), \VH) = 44 1 0) and \HV) = 44 1 0)- We see 
clearly that several processes may lead to the simultaneous firing of the trigger detectors 
(which observe the spatial modes) e and d. The signal photons enter the polarising beam 
splitters. Their action can be described by the following relations 

Xy = cos(6'j)x+ + sin(6'j)xL 
x^'fj = — sin(6'j)x^^ + cos(^^j)x^, (4.5) 

with X = a, b;i = 1,2 respectively and +, — denoting the output spatial modes. 



4.3 Quantum predictions 

The probabilities of various two-particle processes that may occur at the spatially separated 
observation stations, under the condition of both trigger detectors firing simultaneously, are 



on 



given by: 



P(la+, la_; 06+, 06_) = P(2a+, Oa_; 06+, 06_) 
= P(Oa+, 2a_; 06+, 06_) = P(Oa+, Oa_; 16+, 16_) 
= P(Oa+, Oa_; 26+, 06_) = P(Oa+, Oa„; 06+, 26_) = ^, 
P(la+, Oa_; 16+, 06_) = P(Oa+, la_; 06+, 16_) = ^ [sin(^i - 62)]^ , 
P(la+, Oa_; 06+, 16_) = P(Oa+, la_; 16+, 06_) = ^ [003(^1 - ^2)]^ , (4.6) 

where, for example, P(0a+,0a_; 26+,06_) denotes the probabihty of observing two photons 
at the output 6+, and no photons in the other outputs. 

The Bell correlation function for the product of the measurement results on the signals 
at the two sides of the experiment can be redefined in the way proposed in the previous 
chapter, i.e., all standard Bell-type events are assigned their usual values whereas all non- 
standard events are assigned the value of one. I.e., if no photons are registered at one side, 
the local value of the measurement is one, if two photons are registered at one side again the 
local measurement value is one. The latter case includes both the event in which the two 
photons end-up at a single detector, as well as those when two detectors at the local station 
fire. Please note, that the experiment considered is a realization of Bell's idea of "event 
ready detectors" (see e.g. [piauser78|| ). Therefore, non-detection events are operationally 



well defined (as the simultaneous firing of the trigger detectors pre-selects the sub-ensemble 
of time intervals in which one can expect the signal detectors to fire). 

The above value assignment method, as it has been in the previous chapter, works per- 
fectly if one assumes that it is possible to distinguish between single and double photon 
detection at a single detector. Therefore, again it is convenient to introduce the parameter 
a. 

The partial distinguishability blurs the distinction between events (at one side) in which 
there was one photon detected at say the output ±, and events in which two photons entered 
a the detector observing output it, but the detector failed to distinguish this event from a 
single photon count. In such a case the local event is sometimes ascribed by the local observer 
a wrong value namely —1 instead of 1 (if both photons go to the " — " exit of the polariser and 
the devices fail to inform the experimenter that it is a two photon event, this is interpreted as 
a firing due to a single photon and is ascribed a —1 value). Please note, that if one includes 
less than perfect detection efficiency of the detectors this problem is more frequent and more 
involved (we shall not study this aspect here). 



4.4 Conditions to violate local realism 

Under such a value assignment the correlation function reads: 

^a(^^i,^2) = -^cos(2ei - 262) + ^(1 + 2a), (4.7) 

where a is the numerical value of the distinguishability. When we put into the standard 
CHSH inequality this correlation function it violates the standard bound of 2, only if the 
distinguishability satisfies a > ~ 0.948. Such values are definitely beyond the current 
technological limits. As the efficiency of real detectors makes this problem even more acute, 
one has to propose a modification of the experiment that gets rid of this problem. 



in 



4.5 Proposal of modification of the experiment 



Therefore, in front of the idler detector Te we propose to put polarising beam splitter that 
transmit only vertical polarisation whereas in front of the idler detector Td one that transmits 
only horizontal polarisation. This further reduces the relevant terms in our state, i.e. those 
that can induce firing of the trigger detectors, to the following ones: 



1^) = {z{aWv + blb^n) + - <^Wh)) 44|0). (4.8) 

Again we have normalised the above state. 

Using the above formula we can calculate the probabilities of all possible processes in 
this interferometric set-up, conditional on firings of the two trigger detectors: 



P(la+, la_; 06+, 06_) = f [cos(2^i)]^ , 
P(2a+, Oa_; 06+, 06_) = P(Oa+, 2a_; 06+, 06_) = | [sin(2^i)]^ , 

P(Oa+, Oa_; 16+, 16_) = § cos(2e2)^ 
P(Oa+, Oa_; 26+, 06_) = P(Oa+, Oa_; 06+, 26_) = | [sin(2^2)]^ , 
P(la+, Oa_; 16+, 06_) = P(Oa+, la_; 06+, 16_) = ^ [sin(0i - 62 



P(la+, Oa_; 06+, 16„) = P(Oa+, la_; 16+, 06_) = ^ [cos(^i - 62)^ . (4.9) 
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Under the earlier defined value assignment the correlation function for the current version 
of the experiment reads: 

^a(^l,^2) =-|cOs(201-202) 

+ 1(1 - a) [(cos2^i)2 + (cos 2^2)^] + |a. (4.10) 

When such a correlation functions are inserted into the CHSH inequality one has: 

-2 < -| [cos 2(^1 - 62) + cos 2(^1 - 9'^) 
+ cos 2{e[ - 62) - cos 2{e[ - 9'2)] 
+ 1(1 - a) [(cos2^i)2 + (cos2^2)'] + |« < 2. (4.11) 

Please note that some of the terms of the correlation function which depend only on one 
local angle cancel upon insertion into CHSH inequality. 

For a = 1 (perfect distinguishability) the middle expression in ([4.11|) reaches 2.16569, 



i.e. we have a clear violation of the local realistic bound. This maximal violation occurs at 
angles (in radians) 2^i = -1.30278, 2e[ = -2.87435, 2^2 = 1.05326, 26'^ = 2.62386. What 
is more interesting, for a = 0, i.e. for a complete lack of distinguishability between two and 
single photon events at one detector, the expression in ( [4.11| ) reaches a value which is not 
much lower, namely 2.11453. This can be reached for the orientation angles 26i = 0.0837317, 
2e[ = -1.0749, 2^2 = 3.05769, 26'^ = 4.21568). 



Q1 



4.6 Conclusions 



Therefore we conclude that the proposed modification of the entanglement swapping exper- 
iment, despite the unwanted additional events due to the impossibility of controlling the 
spontaneous emissions at the two separate sources, makes it possible to consider it as test of 
Bell inequalities. The standard configuration can serve as a test of local realism only under 
the condition of extremely high distinguishability between two and single photon counts. 

Finally let us mention that the proposed modification in the configuration enables one 
to observe, in the event ready mode, a bosonic interference effect similar to the Hong-Ou- 
Mandel dip | |Hong87| ]. It is described by the first four formulas of ( |4.9| ). E.g. if 6i = tc/A, no 



coincidences between firings of the two detectors of the station a are allowed. All two photon 
events at this station are, under this setting, double counts at a single detector. Thus, we 
have two very interesting non-classical phenomena in one experiment. 
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Part II 

New generalised Bell inequalities and 
GHZ paradoxes for quA^its 



In the first chapter of this part the approach to the Bell theorem employing Bell inequal- 
ities is generahsed to GHZ correlations for which each local observer is allowed to use more 
then two settings of his or her measuring apparatus. 

In the second chapter the functional Bell inequality is derived. Although the functional 
inequalities for two qubits are of less practical importance they seem to be the first step 
towards an analytic search for the critical visibility of two-qubit sinusoidal interference fringes 
which violate local realism. 

In the third chapter we derive the series of GHZ paradoxes for N and N + 1 maximally 
entangled quA'"its observed via unbiased multiport beamsplitters. 



Chapter 5 



Wringing out better Bell inequalities 
for GHZ experiment [3] 



5.1 Introduction 



Greenberger-Horne-Zeilinger correlations | |Greenberger89[ | lead to a strikingly more direct 



refutation of the argument of Einstein Podolsky and Rosen (EPR), on the possibility of 
introducing elements of reality to complete quantum mechanics [[Einstein35|| , than consider- 



ations involving only pairs of qubits. The EPR ideas are based on the observation that for 
some systems quantum mechanics predicts perfect correlations of their properties. However, 
in the case of three or more qubits, in the entangled GHZ state, such correlations cannot be 
consistently used to infer at a distance hidden properties of the qubits. In contradistinction 
to the original two qubit Bell theorem, the idea of EPR, to turn the exact predictions of 
quantum mechanics against the claim of its completeness, breaks down already at the stage 
of defining the elements of reality. 

The reasoning of GHZ involved perfectly correlated qubit systems. However, the actual 
data collected in a real laboratory would reveal less than perfect correlations, and the imper- 
fections of the qubit collection systems would leave many of the potential events undetected. 
Therefore the original GHZ reasoning cannot be ever tested in the laboratory, and one is 
forced to make some modifications (already, e.g., in ||Greenberger9CI|| ) . 

To face these difficulties several M qubit Bell inequalities appeared in the literature 
Mermin9TT5 , Roy9 1 , Ardehali^ , Belinskil?B| , Ziukowski93b |. All these works show that quan- 



tum predictions for GHZ states violate these inequalities by an amount that grows exponen- 
tially with M. The increasing number of qubits, in this case, does not bring us closer to the 
classical realm, but rather makes the discrepancies between the quantum and the classical 
more profound. 

The study of three or more qubit interference effects does not seem to be a good route 
towards a loophole free test of the hypothesis of local hidden variables. However, to interpret 
the results of such experiments]^ one should know the borderline between the quantum and 
the classical (i.e., local realism). According to current literature (with the exception of 
^ukowski93b|| ) we enter the non-classical territory when the fringes in a M qubit interference 
experiment have visibilities higher than 22(^~'^). The principal aim of this chapter is to show 

^ First observation of GHZ correlations has been already reported [|Bouwmeester99[ . 



that, if one allows each of the local observers to have three measurements to choose from 
(instead of the usual two), the actual threshold is lower (for M > 3). 



5.2 Geometrical method of finding Bell inequalities 

Let us first explain the method that will be used in the next two sections (it is called a 
geometrical method) [[Zukowski93b ] on an example of two maximally entangled qubits. The 
method will be presented in the case of stochastic local hidden variables. Application to 
deterministic case is straightforward [[Zukowski93b| 



We consider the state ( p.l|) on which two spatially separated observers a, b measure 
dichotomic observables On{(pn) with eigenvalues ±1 referring to eigenstates 

|±,(/)„)„ = i=(|0)„±e(**")|l)„) (5.1) 

controlled by knobs (local settings) 0„ (n = a,b). 

Let us further assume that in the experiment observer a chooses between, say, Na settings 
of the local apparatus denoted by 0^, i = 1, 2, . . . , Na, and the observer b chooses between 
Nf, settings denoted by 0^ j = 1, 2, . . . , Nb. The quantum prediction for observer a to obtain 
the result m = ±1, and observer b to obtain m' = ±1 is equal to 

P^Um, m'\<Pl, <pi) = \{l- V^mm' cos(0j, + 4)) , (5.2) 

where V2 takes into account the possible less than perfect visibility (see the discussion of the 
notion of visibility in Introduction - equation ( |2.14|) and the discussion below). 

The question is if this set of probabilities is reproducible by local hidden variables, i.e., 

by 

Pw(m,m'|C,0^,) = !^dXp{X)Pa{m\\,<i^^)P,{m'\X,<Pi) (5.3) 

(see also ( p.l3| )). 

The set of quantum probabilities ( ^.21) as well as the set of local hidden variable proba- 
bilities ( |5.3| ) can be treated as the components of 4 x iVa x A*";, dimensional real vectors P^m 
and Phv with components 



PQ\j{mm'; ij) = ^ (l ~ V2mm' cos(0^ + 
PHvimm';ij) = [ t/Ap(A)P,(m|A, 0^)Pfc(m'|A, 0^,) (5.4) 



where i,j,m,m' enumerate the components of the vectors. Both Pq%i and Phv can be 
considered as vectors belonging to a real Hilbert space with the following scalar product 

{F\G) = El=-i Ew=-i El\ Efii Fimm'; tj)G{mm'; ij), (5.5) 

where F and G are arbitrary vectors from this space. 

In every real Hilbert space any two vectors F, G are equal, i.e. F = G, if and only if 

{F\G) = \\F\\^ =\\G\f. (5.6) 
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Thus, if one has two vectors F, G and one knows the norm of, say, the vector F, and 
{F\G) < then F ^ G. This simple observation is especially useful for us because 

we can always calculate the norm of the vector -Pq^, which is not possible for the vector 
Phv- However, as we will see next, one can estimate the scalar product of -Pqm ^^^1 Phv- 

If {Pq\j\Phv) < II-PqmIP one has the Bell inequality. 

As the simplest example let us show how one can obtain within this approach the standard 
threshold value of the visibility (obtained using the CHSH inequality) for the situation in 
which both observers have two local settings to choose from. To this end let us consider the 
case in which Na = Nf, = 2 and let us choose (f)^ = 0, (pl = |, and (pl = — |, (pl = +j. For 
such local settings one has 

P^i^(m,m'; 1,1) = ^(1 -mm'-^) 
PQliim, m- 1, 2) = i(l - mm'-^) 
^(5m("^,"^';2,1) = i(l -mm'-^) 

P^i,(m, m'; 2, 2) = 1(1 + mm'^) (5.7) 

and the norm of | |-Pqm| P equals 1 + \V2 ■ The next step is to write down explicitly the scalar 
product {Phv\Pq\i) 

{PHv\P^li) = J t/Ap(A)[l-^\/2{/a(0i,A)J,(0^,A) + J,(0^,A)J,(02,A) 
+/,(02, A)/,(0i, A) - IM, \)h{4>l A)}], (5.8) 

where, for instance, Ia{4>'i) = Z]m=-i "^-Pa("^|A, 0^). We see that |/a(0^)| < 1 and this 
immediately implies that the modulus of expression in the curly brackets never exceeds 2. 
Thus, 



(p^v|p<yy <i + v^2^. (5.9) 



Comparing ( [5. 9] ) with the norm of ||Pq|//|P one obtains that the necessary condition for Phv 
to be equal to Pqm is ^2 < 



5.3 Quantum and local realistic description of the ge- 
danken experiment 

We consider a source emitting M qubits each of which propagates towards one of M spatially 
separated measuring devices. The generic form of a GHZ M qubit state is 

|v^(M)) = -^(|0)i . . . |0)m + |l)i . . . |1)m). (5.10) 



Let as assume that the operation of each of the measuring apparatus is controlled by a knob 
which sets a parameter (pi, and the l-th apparatus measures a dichotomic observable Oi{(t)i) 
with two eigenvalues ±1 and the eigenstates defined by |±,0i)i = ^ ± e^*'^'''!!);) . The 
quantum prediction for obtaining specific results at the M measurement stations (for the 
idealised, perfect, experiment) reads 



QM 

1 

2M 



l + nfiir,cos (Ei 
[ri equal to —1 or +1). The GHZ correlation function is defined as 



(5.11) 



EW(0i,...,0m) 



n, 



(5.12) 



and in the case of quantum mechanics, i.e. for p(*^) = PqJ, it reads EqJj- = cos{J2i=i4>i)- 

Local realism implies the following structure of probabilities of specific results (compare 
with (|1|)) 



Pffviri, . . .,rM\4>i 

M 



.4>m) 



(5.13) 



where Pj(r;|A,0/) is the probability to obtain result r/ in the l-th apparatus under the con- 
dition that the hidden state is A and the macroscopic variable defining the locally measured 
observable is set to the value 0/0. The locality of this description is guaranteed by indepen- 
dence of Pi on 0j for all i I. 



5.4 Derivation of Bell inequalities via the geometrical 
method 

We shall now derive a series of inequalities for the M qubit GHZ processes based on the 
already mentioned geometric method [[Zukowski93b|| . 



We assume that each of the M spatially separated observers has three measurements to 
choose from. The local phases that they are allowed to set are (f)\ = rr/G, (pl = 7r/2, 03 = 57r/6 
(for the first observer) and for all the other M — 1 observers they are = 0, 02 = vr/3, 03 = 
27r/3, (« = 2,...,M). 

Out of the quantum predictions for the M qubit correlation function at these settings 
one can construct a matrix endowed with M indices 

<i(0,\, • • • 0i!fj = COS {EtU <Pt) = (5.14) 

(ifc = 1,2,3). 



^ We present our results only for the case of deterministic local hidden variable theories, i.e., P;(n|A, 0;) 
1,0, but generalisation to stochastic ones is obvious. 
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All that we know about local hidden variable theories is that their predictions (for the 
same set of settings as above) must have the following form: 



^HV Wii^ ■ ■ ■ 



H 



(M) 



Jl,...,JM ' 



where 



(5.15) 
(5.16) 



4(A,0fJ = E., r,P,(rfc|A,0fJ. 

Of course, in the case of a deterministic local hidden variables theory Jfc(A,0) = ±1. 
is our test matrix. Please note, that all one knows about is its structure. 

The scalar product of two real matrices is defined by 



(5.17) 



Our aim is to show the incompatibility of the local hidden variable description with the 
quantum prediction. To this end, we shall show that, for two or more qubits, 

(g(M)^^(A/)) 

<2*^-iv^< ||g(*^)||2 = . (5.18) 
First, we show that HQ'-^^^IP = 3^^/2. This can be reached in the following way: 

)||2 



= + exp (2EEi0O1} 

--3'^/2 + Re (nfciiE?,=iexp(2^0: 



k 

ik' j ' 



(5.19) 



where Re denotes the real part. Since I]f=i e**'' i)(2/3)7r _ ^^iq last term vanishes. 

The scalar product {H^^'^\Q^^^^) is bounded from above by the maximal possible value 

of 



S 



(M) 



E 

i\,...,iM 



M 



M 



coME< n^KA,<) 



\k=l 



1=1 



and for M > 2 



S 



(M) 



< 2 



M~l 



V3. 



To show ( ^.21|) , let us first notice that 



S 



(M) 



Re 



M 3 



n E4(A|0t)exp«,)) 



k=l ifc=l 

M, one has e^'^i = e*'*^'"^^/^''" whereas for k = 1, e^'^i 



For A; = 2, . 
since /(A|-) = ±1, the possible values for 



J2 h{X\<j>l) exp{z<Pl: 



(5.20) 
(5.21) 

(5.22) 

gi{7r/6)gi[{«-l)/3]7r^ ThuS, 

(5.23) 



h=i 



An 



are 0, ±2e^''/^ ±2e"^(''/^), or finally ±2e^(''/^), whereas for = 2, . . . , M the possible values 
of 



^4(A|0^Jexp(20, 



(5.24) 



2 = 1 



have their complex phases shifted by 7r/6 with respect to the previous set; i.e., they are 0, 
_^2e^(2V3)^ ±2, or finally ±2e''^'^'^\ Since Iz^Wf^^zll < 2*^ and the minimal possible overall 
complex phase (modulo 2tt) of z^Yl^^^2^k is vr/6, one has Re{z^ Yl^^2 ^k) — 2*^ cos(7r/6). 
Thus inequalities ( 5.21|) and ( ^.18 ) hold. 

The left inequality of ( |5.18| ) is a Bell inequality for the M qubit experiment. If one 
replaces H^^^ by the quantum prediction Q^^'^^ (compare ( ^.14D ) the inequality is violated 
since 

(gW,gW) = 2_>2*^-iv^, (5.25) 

i.e., (pl8D is violated by the factor (3/2)*7V3 (compare ||Mermin90b| |). 



5.4.1 Critical visibility and quantum efficiency of detectors 

The magnitude of violation of a Bell inequality is not a parameter which can objectively 
define to what extent local realism is violated. It is rather the visibility of the M qubit 
interference fringes which can be directly observed. Further, the significance of all Bell- 
type experiments depends on the efficiency of the collection of the qubits. Below a certain 
threshold value for this parameter experiments cannot be considered as tests of local realism. 
They may confirm the quantum predictions but are not falsifications of the hypothesis of 
local hidden variables. Therefore we will search for the critical minimal visibility of M qubit 
fringes and collection efficiency, which do not allow anymore a local realistic model. 

In a real experiment (under the assumption that quantum mechanics gives idealised, but 
correct predictions), the visibility of the M qubit fringes, V2(M), would certainly be less than 
1. Also the probability of registering all potential events rj2{M) would be reduced by the 
overall collection efficiency. If one assumes that all M local apparata have the same collection 
efficiency rj, and takes into account that these operate independently of each otherQ, one can 
model the expected experimental results by 

nl^5(ri,...,rM|0\...,0^) 
= 77^(1) ''(l + r2(M)nfiincosE^ii0'=) . (5.26) 

The full set of events at a given measuring station consists now of the results +1 and —1, 
when we succeed to measure the dichotomic observable, and a non-detection event (which 
is, in principle observable, if one uses event-ready state preparation [ Yurke9!^ ) for which one 



can introduce the value 0. The local realistic description requires that the probabilities of 
the possible events should be given by 

Pfp^,(mi,...,mM|0\...,0^O 

= Sd\p{\)Y{lUP,{m,\X,<P^) , (5.27) 



■^The parameter rj describes here the efRciency of a single detector. The assuraption of the independency 
of detectors gives 772 (M) = 77*^. 



with rrii = +1,-1 or 0. The local hidden variable correlation function for the experimental 
results (at the chosen settings) is now given by 



M 



k=l 



with 

l'ki\<Pt)= E mkPrim.W^l). (5.29) 

mfe=-l,0,+l 

For deterministic models one has now /^(A, (p^J = 1, 0, —1. 

One can impose several symmetries on P^p^. These symmetries are satisfied by the 
quantum prediction ( |5.26|) , and we can expect them to be satisfied in real experiments, 
within experimental error. The one that we impose here is that: 

For all sets of results, {mi, . . . , mAf}, that have equal number of zeros (one zero 
or more) the probability P^^{mi, . . . , tum) has the same value, and this value is 
independent of the settings of the local parameters {^J^, . . . , 

One can define a function /^/(m) which for m = +1, —1, has the following values: /(±1) = 
±1, /(O) = —1 (compare ||Garg85|| ) and introduce auxiliary correlation function 

Eil,...,iN — S d\p{\) Z]mi,...,m]v=-1,0,+1 

X nf=i[/(mfc)P,(m,|A,0^J] = H\^;y^ . (5.30) 

Since, due to the symmetry conditions, one has, e.g., J2m2=i-i fi''^2)Pi^pt{0, ^2, . . . , = 
0, the following relation results: 

E,,„„,,, = Egl^^_^^^ + [fiO)rm • • • , 0), (5.31) 

where P(0, . . . , 0) is the probability that all detectors would fail to register qubits, and under 
our assumptions it is independent of the settings, and equals (1 — r^)*^. 

The auxiliary correlation function must satisfy the original inequality ( p. 181 ); i.e., one has 

(QW,^W)< 2*^-^73. (5.32) 

However, this implies that 

_2M-iV3-/(0)^^P(0,...,0)g(M) 

< 2^^-'VS - /(0)*^P(0, . . . , O)g(M), (5.33) 

where 



<liM)= E CUr (5-34) 



ll,...,tM 

Therefore, since if x is a possible value for {Q^^^\ Eg]^^) then so is —X, one has 



expt/ 

<2^-iv^-P(0,...,0)|g(M)|. (5.35) 



AO 



Thus, we have obtained Bell inequalities of a form which is more suitable for the analysis of 
the experimental data. 

The prediction ( ^.26 ) leads to the following correlation function 



E. 



QM 
expt 



V 



V2{M)E^^, 



(5.36) 



which, when put into ( |5.35| ) in the place of -E^p^, gives the following relation between the 
threshold visibility, V^*''(M), and the threshold collection efficiency, r]^^, for the M-qubit 
experiment: 

^trM6_^,r^M^ = 2'^^-'VS-\q^M)\il-v'T- (5-37) 



The value of the expression q(^M) can be found in the following way: 

QiM) = En,...,^M COS (Efcll 



Re (En,.. 
= Re {Uti E.,exp(z0, 



Re 



2^iexp (z(M-l)f) 



-2^' sin (M -1 



(5.38) 



5.5 Results 



The threshold value of the visibility of the multi-qubit fringes decreases now faster than in 
the earlier approaches ||Mermin90b|| . For perfect collection efficiency, [i] = 1), it has the 
lowest value, which is 

Vr{M) = V3[-) , (5.39) 

and, if M > 4, it is lower than (:^)^^~^- The specific values for several qubits are V2^{2) = 
76.9%, ^'"(S) = 51.3%, V^'iA) = 34.2%, 1/2*'(5) = 22.8% and Vi'iW) = 3% (see also 
figure (13)), whereas the standard methods lead to V°^'^{2) = 70.7%, ^"''^(3) = 50.0%, 



^2 



(4) = 35.4%, 1/2°''^(5) = 25.0% and ^'"'^(lO) = 4.4%. 



The threshold efficiency of the qubit collection also decreases with growing M (see figure 
I)), and for perfect visibilities it reads r/*''(2) = 87.0%, r/*^(3) = 79.8%, t]^''{A) = 76.5%, 



r7*'"(5) = 74.4% (here the number in the brackets indicates the number of entangled particles 



M). The gain over the inequalities | [[Vlermin90b| | is in this respect very small, and begins 



again at M = 4. However, for very big M the critical efficiency is close to ^ (compared with 



^ for [|Mermin90b|] ). 



5.6 Conclusions 

We conclude that for the original GHZ problem (four qubits) one should rather aim at making 
experiments which allow for three settings at each local observation station. Surprisingly, the 
measurements should not be performed for the values for which we have perfect GHZ-EPR 
correlations (i.e the values for which the correlation function equals to ±1). 
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Figure 5.1: The values of critical visibility V2^{M) versus the number of qubits M. 
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Figure 5.2: The values of critical efficiency rf'^{M) versus the number of qubits M. 
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Chapter 6 



Bell inequality for all possible local 
settings [4] 



6.1 Introduction 



The Bell theorem is usually formulated with the help of the Clauser-Horne [ UlauserT^ or the 
CHSH inequality [|Clauser69|| . These inequalities are satisfied by any local realistic theory 
and are violated by quantum mechanical predictions. They involve two apparatus settings at 
each of the two sides of the experiment. However, generalisation to more than two settings 



at each side are possible ||Garuccio80|, |Braunstein89| , |Zukowski93b| , |Gisin99 . 

There are several motivations for such generalisations. First of all, new Bell inequalities 
may be more appropriate in some experimental situations, e.g., the chained Bell inequalities 
can reveal violation of local realism for the Franson type experiment |[Aerts99|] . Also, the 



academic question, why only two settings at each side, is that always necessary, is interesting 
in itself. Further, many of the currently performed quantum interferometric Bell tests did 
not involve stabilisation of the interferometers at specified settings optimal for the standard 
Bell inequalities, but rather involved sample scans of the entire interferometric patterns. 
Thus it is useful to have inequalities that are directly applicable to such data. 

Here we present a Bell-type inequality that involves all possible settings of the local 
measuring apparatus for a pair of two qubits, which is always equivalent to two spin ^ 
particles. The method applied is a development of the one given in [|Zukowski93b|| . However, 
here we do not restrict ourselves to pairs of coplanar settings (in the meaning appropriate 
for two Stern-Gerlach apparatuses). 

Our method has two characteristic traits. The first one is that it indeed involves the entire 
range of the measurement parameters. By this, e.g., it distinguishes itself from the limits 
of infinitely many settings at each side of the so-called chained inequalities ||Braunstein89 



in which not every pair of possible settings is utilised. The second one is that the method 
involves the quantum prediction from the very beginning. As we shall see the quantum 
prediction determines the structure of our Bell inequality. 



6.2 Quantum mechanical and local realistic description 
of the gedanken experiment 



As usual one has a source emitting two qubits each of which propagates towards one of two 
spatially separated observers a and b. The qubits are described by the maximally entangled 
state (^. 

Let as assume that every observer has a Stern-Gerlach apparatus, which measures the 
observable n ■ a, where n = a,b, n is a unit vector representing direction at which observer n 
makes a measurement and a is a vector the components of which are standard Pauli matrices. 
The family of observables n ■ a covers all possible dichotomic observables for a two qubit 
system, endowed with a spectrum consisting of ±1. 

In each run of the experiment every observer obtains one of the two possible results of 
measurement, ±1. The probability of obtaining by the observer a the result m = ±1, when 
measuring the observable a ■ a, and the result m' = ±1 by the observer b, when measuring 
the observable a ■ b is equal to 

PQ%j{m,m'\a,b) = j{l — V2mm'a- b). (6.1) 

where < V2 < 1 stands for the visibility. 

The structure of local hidden variables gives 

-Pffy("^, m'|a, 6) = (iAp(A)Pa("^|A, a)Pb(m'|A, 6), (6.2) 



with the standard meaning of the used symbols (see (p.l3|)). 



6.3 Derivation of the inequality via the geometrical 
method 

To apply the geometrical method we must define appropriate Hilbert space. Because we deal 
with functions PQ|^(m, m'|^a, 0a, 6*6, 4>b) and PHvijri) fn'\9a, 4>a, Ob, 4>b) that depend on discrete 
numbers m, m' and continuous variables 9n, (pn, where n = (sin 9n cos sin 6'„ sin 0„, cos 9n) 
it is convenient to define the scalar product of certain two real functions / and g as 

{f\9) = E E 

m=—l m'=— 1 

J dria J dni,f{m, m'; 9a, (t)a, 9b, (j)b)g{m, m'; 9a, (pa, 9b, (pb), (6.3) 

where dfln = sin9nd9nd(f)n is the rotationally invariant measure on the sphere of radius one. 
Our known vector is Pqm, whereas the test one is Phv (compare with the section describing 
geometrical method). 
One has 

I I pV2 ||2_/pV2 I pVi \ 
I I Q Af 1 1 — VQM\-^QMI 

= (2^)^ + Vi^. (6.4) 

A7 



To estimate the scalar product {Pqm\Phv) one has to use the specific structure of proba- 
bihties that are described by local hidden variables (LHV) (|6.2|). Since Phv is a weighted 
average over the hidden parameters one can make the following estimate 



{Pqm\Phv) < max 



AeA 



J2 J dila J dVtbPa{mWa) 



m,m'=— 1 



xPfe(m'|A, — mm'V2d ■ b) 



Therefore ( |6.8| ) reads 



^ k=l-' 



(6.5) 



Since 

YL=-i Pa{m\\ a) = Eln'=-i Pb{m'\X, b) = 1, (6.6) 
the first term of ( |6.5| ) satisfies 

^ 51 / d^a J dQbPa{m\X,a)Pb{m'\X,b) 

m,m'=— 1 

= (2vr)^ (6.7) 
We transform the other term of (B.SI) to a more convenient form 



1 ^ r r 

- ^ y dila J d^lbrnm'Pa{rn\X,a)Pb{rn'\X,b)Va-b 

m,m'=—l 

\j d^aj dnjaid, X)h{b, X)V2d ■ b, {6.1 



where 

In{n, X) = Em=-i rnPn{m\X, n), (6.9) 

and one has |/„(n, X)\ < 1 {n = a, b). 

The scalar product of two three dimensional vectors a and b that appears in (|6.8| ) can be 
written as a ■ 6 = Y%=i ak{Oa, 4>a)bk{0b, (pb), where 

n = (ni,n2,n3) 

= (sin 6'n cos0ri, sin 6'n sin0„, cos6'„). (6.10) 



X J d^ibib{ebAbA)hk{ebAb). (6.11) 

We notice here that our expression is a sum of three terms, each of which is a product of 
two integrals. 

The functions in ( |6.11D are square integrable, i.e. integrals 



j dl]„|/„(^„,0„,A)p (6.12) 



/IS 



and 

diln\nk{0n,4>n)\'^ (6.13) 



/■ 



exist (we remind that \In{On, 'Pn, A)| < 1 which guarantees the existence of the first integral). 
This allows us to use formalism of Hilbert space of square integrable functions on the unit 
sphere, which we denote as L'^{S^). 

The functions nk{9n,4>n) fulfil the orthogonality relation J dflnnk{On, (pnjnii^On, 4>n) = 
^Ski- Thus, if we normalise (i.e. we divide them by their norm, which is y^) we 

can interpret the integral ^^(A) = -y/^ / dQnIn{9n, 4>n, ^)nk{0 k-th coefficient of 

the projection of In{On, 4>n, ^) into a three dimensional subspace of L^(5'^) spanned by the 
(normalised) basis functions \J^nk{0n,4'n) {k = 1,2,3). For later reference we will call this 
space S(3). Therefore (|6.8|) transforms into 

V2^j:at{X)al{X). (6.14) 
fc=i 

Denoting the projection of In{On,(pn,X) into S(3) by /Jl(^^„, 0„, A) and using the Schwartz 
inequality we arrive at 

^E«^(A)«^(A)<^||/^(-,A)||||4l(-,A)||. (6.15) 
k=i 

Therefore, our last step is to calculate the maximal possible value of the norm ||/J|(-, A)||. 
Since the length (norm) of a projection of a vector into a certain subspace is equal to the 
maximal value of its scalar product with any normalised vector belonging to this subspace, 
the norm ||jjl(-, A)|| is given by 



A) 1 1 =max\s\=i[^f^ J dnjri{On,(l)n,X)T,k=iCknk{On,(l)n)], (6.16) 

where c= (01,02,03) and |c|^ = J2k=i(^k — ^- Because [/^(a, A)| < 1 one has 

||jll(-,A)|| <max|,-i=i[yj/rffi„|ELiCfcnfc(^„,0O|]. (6.17) 

Every vector c can be obtained by a certain rotation of the versor z. Such a rotation is 
represented by an orthogonal matrix O belonging to the rotation group S0{3). Therefore, 
( |6.17D can be rewritten as 



\\l\\{;X)\\<maXoyiJdnn\Oz-n{en,(l)n)\], (6.18) 

where the maximum is taken over all possible rotation matrices O. Since \0z ■ n(6'„,0„)| is 
the modulus of the scalar product of two ordinary three dimensional vectors, it is equal to 
\z ■ 0~^n{9n,(l>n)\- An active rotation of the vector n is equivalent to a (passive) change of 
the spherical coordinates. Utilising the fact that the measure dQn is rotationally invariant 
we see that 

||/iJ||< y'dfi„|J|-cos^„| = 27rJ|-. (6.19) 



An 



Therefore (|6.15| ) is not greater then i(27r)^, which with 
following inequalities 

WPqIiW = (2^)' + > (2vr)2 + ^(27r)^ > {P^UPhv). (6.20) 

6.4 Results 

The inequality ( |6.2U| ) is violated by quantum predictions provided that the visibility V2 is 
higher then 75%. Please notice that the right hand inequality is a form of a "functional" 
Bell inequality. It simply gives the upper bound for the value of a certain functional defined 
on the local realistic probability functions Phv- The left hand inequality shows that the 
insertion of Pqm i^i^o the functional Bell inequality leads to its violation provided V2 > 0.75. 
The characteristic trait of our functional Bell inequality is that its form is defined by the 
quantum prediction Pq^. 



6.4) and (|6.7|) gives us the 



6.5 Conclusions 



The threshold visibility for two qubit interference to violate the inequality (|6.20|) is lower 



than in the case of coplanar settings | Zukowski93b |, for which the critical visibility is 
Also, it is lower than the one given recently by Gisin ||Gisin99|| . For his inequalities involving 
arbitrary many settings the threshold visibility equals V2 = f . The chained inequalities 
paruccio80|| , |Praunstein89|] , for evenly spaced settings, with the number of settings going 
to infinity, have the property that the critical visibility approaches 1 in the limit of infinitely 
many settings. 

Taking into account the fact that the necessary condition for the existence of a local 
hidden variable model for two local settings at each side of the experiment implied by the 
CHSH inequality is V2 < the critical visibility of any hidden variable model that aims 
at reproducing quantum correlations for the gedanken experiment described above for any 
number of local settings (infinite or not) cannot be greater than Therefore, the numerical 
value of the critical visibility obtained by the above inequality is overestimated. This is due 
to the fact that this inequality is only a necessary condition for the existence of local hidden 
variables. However, the inequality ( |6.20[ ) may be the first step towards finding the threshold 
visibility of the local hidden variable model reproducing quantum mechanical predictions 
for all positions of measuring apparatus. The presented inequality also solves the academic 
problem of finding a Bell inequality that involves all possible settings of the local apparata. 

In ||Zukowski93b | Zukowski has shown that for GHZ states involving four or more par- 
ticles (and employing all coplanar settings) the functional Bell inequality approach leads to 
much more stringent conditions on the critical visibilities than the approaches of Mermin 



Mermin90b|| and Ardehali [|Ardehali92|| . However, the extension of the presented approach 
to GHZ states does not give any improvement with respect to the standard one. This is due 
to the fact that original GHZ paradox is obtainable for co-planar settings only. 



Chapter 7 



Greenberger-Horne-Zeilinger 
paradoxes for quA^its [5,6] 



7.1 Introduction 

The Greenberger-Horne-Zeilinger correlations, discovered in 1989, started a new chapter in 
the research related with entanglement. To a great extent this discovery was responsible for 
the sudden renewal in the interest in this field, both in theory and experiment. All these 
developments finally led to the first actual observation of three qubit GHZ correlations in 
1999 I Bouwmeester99 1 , and as a by product, since the experimental techniques involved were 



of the same kind, to the famous teleportation experiment ||Bouwmeester97| 



In this chapter we would like to examine whether GHZ-type paradoxes exist also in the 
case of correlations expected in gedanken experiments involving multiport beam splitters 
[[Klyshko88| , |Zeilinger9^ , |Zeilinger94|| , i.e. for a specific case of non dichotomic observables 
(which have properties distinctive to the dichotomic ones ||GleasonFr| , |Bell66| , ^ocEen^). 
To this end, we shall study a GHZ-Bell type experiment in which one has a source emitting 
M quA^its in a specific entangled state of the property, that the quA^its propagate towards 
one of M spatially separated non conventional measuring devices operated by independent 
observers. Each of the devices consists of an unbiased symmetric multiport beam splitter 
[ [Zukowski97b|| (with A^ input and A^ exit ports), A^ phase shifters operated by the observers 
(one in front of each input), and A^ detectors (one behind each exit port). 



7.2 Unbiased multiport beamsplitters 



An unbiased symmetric 2A^-port beam splitter is defined as an A^-input and A^-output inter- 
ferometric device which has the property that a beam of light entering via single port is evenly 
split between all output ports (see (|7.1|) ). I.e., the unitary matrix defining such a device has 
the property that the modulus of all its elements equals An extended introduction to 

the physics and theory of such devices is given in ||Zukowski97b|| , and therefore the reader not 
acknowledged with those concepts is kindly asked to consult this reference. Multiport beam 
splitters were introduced into the literature on the EPR paradox in ||Klyshko88| , |Zeilinger93| , 
Zeilinger"94 ] in order to extend two qubit Bell-phenomena to observables described as opera- 
tors in Hilbert spaces of dimension higher than two. In contradistinction to the higher than 
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Figure 7.1: The schematic picture of an 2N port. There are N phases in front of each input 
port and detectors behind each output ports. 



1/2 spin generahsations of the Bell-phenomena |[Mermin80| , |Garg82| , [MerminS^ , |Ardehali91| , 
AgarwalQ^ , |W6dkiewicz94| , |W6dkiewicz95| ], this type of experimental devices generalise the 
ideas of beam-entanglement |Pornc84 Zukowski88 , Horne89 , [Rarity9"0[| . Unbiased symmetric 
multiport beam splitters are performing unitary transformations between "mutually unbi- 
ased" bases in the Hilbert space Pchwinger60| , [lvanovic8T| , [Wooters86| . They were tested in 



several recent experiments ||Mattle95| , [Reck96|] , and also various aspects of such devices were 
analysed theoretically [[Reck94| , |Jex95|| . 

We shall use here only multiport beam splitters which have the property that the elements 
of the unitary transformation which describes their action are given by 



N 

m,m' 



1 

IN ) 



(7.1) 



where = exp(z-^) and the indices m, m' denote the input and exit ports. Such devices 
were called in ||Zukowski97b|| the Bell multiports. 



7.3 Quantum mechanical predictions 

One assumes that the initial M quA^it state that feeds M spatially separated multiports, 
each of which has A^ inputs and A^ outputs, has the following form: 

1 N M 

= ^ E n 1^)^ (7-2) 

where |m)/ describes the /-th quA^it being in the m-th beam, which leads to the m-th input 
of the Z-th multiport. Please note, that only one quA^it enters each multiport. However, each 
of the quA^its itself is in a mixed state (with equal weights), which gives it equal probability 
to enter the local multiport via any of the input ports. 

The state ( [7.2| ) seems to be the most straightforward generalisation of the GHZ states 
to the new type of observables. In the original GHZ states the number of their components 
(i.e., two) is equal to the dimension of the Hilbert space describing the relevant (dichotomic) 



degrees of freedom of each of the quA^its. This property is shared with the EPR-type states 
proposed in ||Zukowski97b|| for a two-mult iport Bell-type experiment - in this case the number 
of components equals the number of input ports of each of the multiport beam splitters. We 
shall not discuss here the possible methods to generate such states. However, we briefly 
mention that the recently tested entanglement swapping ||Zukowski93a| , |Zukowski95| , |Pan98|| 
technique could be used for this purpose. 

As it was mentioned earlier, in front of every input of each multiport beam splitter one 
has a tunable phase shifter. The initial state is transformed by the phase shifters into 



N M 



IV^(M)') = ^ E n«^P(^C)l^)z, (7.3) 



N m=l 1=1 

where 0™ stands for the setting of the phase shifter in front of the m-th port of the l-th 
multiport. 

The quantum prediction for probability to register the first photon in the output ki of 
an 2N - port device, the second one in the output ^2 of the second such device and the 
M-th one in the output of the M-th device is given by: 



PqAiiki, . . . , kM\(pi, 



1\M+1\X^N Y^M Am\j-fM .,(m-l){fcn-l) 12 



E^=iexp(zErii0r)nr=i7]v' 



'J_\A/+1 



;iV + 2EL„.cos(EfiiA$-,-')], (7.4) 

where A$J"^™ = 0™ — 0™' + ^{ki — l){m — m'). The shorthand symbol 0^ stands for the 
full set of phase settings in front of the k-th multiport, i.e. 0^, 0|, . . . , 0^. 

7.3.1 Bell number assignment 

To efficiently describe the local detection events let us employ a specific value assignment 
method (called Bell number assignment; for a detailed explanation see again |[Zukowski97b|| ), 
which ascribes to the detection event behind the m-th output of a multiport the value 7]v~^, 
where 77V = exp(i^). With such a value assignment to the detection events, the Bell-tj^e 
correlation function, which is the average of the product of the expected results, is defined 
as 



E(0i,---,0A/) = 

= J:t-,k,,=i nfii l'^-'P{ki, ■ ■ ■ , A;m|0";, ■ ■ ■ , 0m) (7.5) 

and as we shall see for the quantum case it acquires particularly simple and universal form^. 

The easiest way to compute the correlation function for the quantum prediction employs 
the mid formula of ( [7.4|) : 

^QAf(01,---,0Af) = 
^This is the main purpose for using this non-conventional value assignment. 



KA 



>^ lll=llN — 



= iir^' E^,^,=i exp Eti(C - 

xnfiiE£=i7S'-'^^"'+'^ (7.6) 

Now, one notices that Y^k=i In'' differs from zero (and equals to N) only if m — 

m' + 1 = 0, modulo N. Therefore we can finally write: 



Eqm{4>i, - ■ ■ Am] 

N ^ 



where (j)^'^'^'^ = — 0™"*"^ and the above sum is understood modulo N, which means that 

One can notice here a striking simplicity and symmetry of this quantum correlation 
function (|7.71 ). It is valid for all possible values of M (number of quiVits) and for all possible 



values of > 2 (number of ports). For = 2, it reduces itself to the usual two qubit, 
and for N = 2, M > 2 the standard GHZ type multi-qubit correlation function for beam- 
entanglement experiments, namely cos(X]fii <p]''^) | |Mermin90b| ]. The Bell - EPR phenomena 
discussed in ||Zukowski97b|| are described by (|7l7| ) for M = 2, > 3. 



Even for = 2, M = 1 the function ( [7.7D describes the following process. Assume that 
a traditional four-port 50-50 beam splitter, is fed a single photon input in a state in which 
is an equal superposition of being in each of the two input ports. The value of (|7.7| ) is the 
average of expected photo counts behind the exit ports (provided the chck at one of the 
detectors is described as +1 and at the other one as —1), and of course it depends on the 
relative phase shifts in front of the beam splitter. In other words, this situation describes a 
Mach-Zehnder interferometer with a single photon input at a chosen input port. For A^ = 3, 
M = 1 the same interpretation applies to the case of a generalised three input, three output 
Mach-Zehnder interferometer described in [|Weihs96|| , provided one ascribes to firings of the 
three detectors respectively 73 = a = exp(«^), and a^. 

7.3.2 Perfect correlations 

The described set of gedanken experiments is rich in EPR-GHZ correlations (for M > 2). 
To reveal the above, let us first analyse the conditions (i.e. settings) for such correlations. 



As the correlation function (|7.7| ) is an average of complex numbers of unit modulus, one has 
\Eqm{(Pi, ■ ■ ■ ,0m) I < 1. The equality signals a perfect EPR-GHZ correlation. It is easy to 
notice that this may happen only if 

M M M 

exp(i <PI'^) = exp(« '^P) = ■ ■ ■ = exp(i ^ (pf'^) = 7^, 
1=1 1=1 1=1 

where k is an arbitrary natural number. Under this condition -E(0i, ■ ■ ■ ,4'm) = In- This 
means that only those sets of M spatially separated detectors may fire, which are ascribed 
such Bell numbers which have the property that their product is 7^. Knowing, which 
detectors fired in the set of M — 1 observation stations, one can predict with certainty which 
detector would fire at the sole observation station not in the set. 



7.4 Paradoxes for + 1 maximally entangled quA^its 



7.4.1 Four tritters 

We shall now present the simplest GHZ-type paradox for such systems. We take N = 3 and 
M = 4. That is we consider now, the experimental situation in which one has the source 
producing the ensemble of four three-state particles (qutrits) described by the state \ip{4:)) 
(compare, (|7.2|) ) that feeds four three-port beam splitters (i.e., tritters [|Zukowski97b| ]). In 
this case the quantum correlation function has the form: 



The (deterministic) local hidden variables correlation function for this type of experiment 
must have the following structure ||Bell64||: 



£^W(01,02,03,04) = / l[hi(pk,\)pWdX. (7.9) 

•^^ fc=l 

The hidden variable function Ik{(f)k, A), which determines the firing of the detectors behind 
the k-th multiport, depends only upon the local set of phases, and takes one of the three 
possible values a, a^, = 1 (these values indicate which of the detectors is to fire), and 
p(A) is the distribution of hidden variables. 

Consider four gedanken experiments. In the first one our observers, each of whom op- 
erates one of the spatial separated devices, choose the following phases in front of their 
three-port beam splitters: 



0i = (0i,0?,0?) = (O,^ 



9 ' 9 



/'i,0i0!) = (o,o,o) = 0'. (7.10) 



In the second experiment, the third observer sets 03 = 0' whereas the other ones set 0. 
We repeat this swapping of the settings procedure in the next two experiments until the 
first observer sets 0' and the other set 0. Quantum mechanics predicts that in all four such 
experiments the correlation function is equal to (i.e. we have perfect GHZ correlations). 
Namely we have 

EqaM, 0, 0, 0') = Eqm{$, 0, 0', 0) 

= Eqm{1 0', 0, 0) = Eqm{$\ 0, 0, 0) = a\ (7.11) 

However, this immediately implies that for any local hidden variables theory that aims 
at describing these phenomena one must have for every A 



4(0',A)nf=w./K0,A) = a2, (7.12) 



and this must hold for all k = 1,2, 3, 4. But, since /z(0, A)^ = a^^ = 1 (where, k represents a 
certain integer), then after multiplying these four equations side by side, one has for every A 



nf=i/K0',A) = a^ (7.13) 

Therefore, if the local hidden variable theory is to agree with the earlier mentioned quantum 
predictions ( [7.11| ), then one must have 



E(0', (/)',(/)', 00 = a - a*. (7.14) 

However, the quantum prediction is EQM{(p',(p',4>'y4>') = 1- Thus we have a GHZ-type 
contradiction that 1 = a*. I.e., hidden variables predict a different type perfect EPR-GHZ 
correlation. I.e. we have a realisation of the GHZ paradox for non-dichotomic observables. 

7.4.2 General case: + 1 maximally entangled quA^its 

We will extend the reasoning to the case when one has M = N + 1 quA^its (described by 
the state of the form ( [7. 2[ ) ) beamed into multiport beam splitters with N input and output 
ports. The quantum prediction for the Bell correlation function is given by (|7.7|) , with the 
appropriate value of M. 

The local hidden variables correlation function must have the following structure 



/nf=i4(V4,A)p(A)rfA, (7.15) 

where k now extends from 1 to M and ipk stands for the full set of settings in front of the 
fc-th multiport, i.e. ipl, ipl, ■ ■ , ip^^'^, and Ikii'k, A) is a hidden variable function depending 
on the local phase settings, which has the property that its value, which can be any integer 
power of 7m-i, indicated which local detector is to fire. 

Now, as it was in the previous case, we must choose appropriate phases for each of 
observers that will be taken in the first experiment. The appropriate choice is the following 
one: 



V'l = • • • = V^M-i = (0, 5,2S,---,{M- 2)S) = ^ 

^M = {0,---,0)=i^', (7.16) 

where 6 = jj^zff- the next M — 1 experiments one applies previously described swapping 

of the settings procedure until the first observer sets ip' and the other ones set ip. For such 
choice of phases the quantum correlation function for every of the M experiment, is equal 
to 7mIi = 7m-i- we have perfect GHZ correlations of the same type for each of 

the experiments). But this implies that, for any local hidden variables theory that aims at 
describing these phenomena one must have, for every A, 

4(VS'',A)nz^fc//(^,A)=7M-i, (7.17) 
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and that this must hold for all /c = 1, • • • , M. However, after multiplying these M equations 
one has: 



UfiiWA) = YM-i, (7.18) 

where we have used the property of the Bell numbers generated by 7m- i, that each of them 
to the M — 1-th power gives 1, and therefore that Ik{ijj,X)^ — 1. Thus, the local hidden 
variable implies that 

EQM{if', - ■ ■ =7m-i- 

However, the quantum prediction is 1. Thus we have the GHZ contradiction that 1 = 7m-i- 
I.e. hidden variables predict a different type perfect EPR - GHZ correlations. In other words 
the EPR idea of elements of reality makes no sense for the discussed experiments, and this 
hold for an arbitrary number of quA^its M, and for suitably related (M — 1), but in principle 
arbitrarily high number of input and exit ports of symmetric multiport beam splitters. 



7.5 Paradoxes for A'' maximally entangled quA^its. 

In this section we show that the above reasoning can be as well applied to the case where the 
number of observers M equals the number of input ports N of the 2A'" port Bell multiports. 

7.5.1 Three tritters 

As the simplest example let us consider the gedanken experiment with three observers each 
of which having tritters (3 input and 3 output ports) as a measuring device. The correlation 
function for such an experiment reads 

-Eqm(01>2,03) 

-5 k=i \ 1=1 ) 

In the first run of the experiment we allow the observers to choose the following settings 
of the measuring apparatus 

7* , TT 27r, 7* 7* 

03 = (O,O,O) = 0' (7.20) 

in the second run they choose 




(7.21) 

whereas in the third run they fix the local settings of their tritters on 

01 = = 03 

02 = 0'. (7.22) 



Now, let us calculate the numerical values of the correlation function for each experimen- 
tal situation. We easily find that for all three experiments this value, due to the special form 
of the correlation function, is the same and reads 

= exp(-zy ) = a\ (7.23) 

i.e., we observe perfect correlations. 

Proceeding in exactly the same way as in the previous section we can write the equations 
that must be fulfilled by local hidden variables for every A, i.e., 

/i(0,A)/2(0,A)/3(0',A) = a2 
/i(0',A)/2(0,A)/3(0,A) = a2 

Ji(0,A)/2(0',A)/3(0,A) = a2 (724) 
After multiplication of the above equations one arrives at: 

3 3 

n 4(0', A) n 4(0, \f = {c^f = 1, (7.25) 

A;=l fc=l 

which can be also written in the following form 

3 3 

n4(0',A) = n4(0,A), (7.26) 

k=l k=l 

where we have used the property of the hidden variable functions, namely that 4(0, A)^ = 
4(0, A)* for every A. However, because i?QM (0', 0', 0') = 1 we must also have 

n4(0',A) = l, (7.27) 

k=l 

which, because of ( |7.26| ), gives 

n4(0,A) = l (7.28) 

for every A. This in turn implies that 

£;qm(0,0,0) = 1, (7.29) 

which means that local hidden variables predict perfect correlations for the experiment when 
all observers set their local settings at cf). However, the true quantum prediction is that 

EQM(0,0,0) = -i. (7.30) 

Therefore we have the contradiction: 1 = — i. 

This contradiction is of the different type than the one derived in the previous section 
although it has been obtained in the similar way, i.e., the perfect correlations have been used 
to derive it- equations ( [7.24|) and (|7. 25|) . Here local hidden variables imply a certain perfect 
correlation, which is not predicted by quantum mechanics. 



7.5.2 General case: N maximally entangled quA^its. 

Now, we employ the above procedure for the case when we have an arbitrary odd number 
of multiports and quiVits N = 2m + 1. As we have seen in the previous section the crucial 
point is to find the proper phases for the multiports. Let us choose for the first gedanken 
experiment the following ones 



= (0, 



TT 



2n 



2rmr 



2m + 1' 2m+ 2m + 1' 



^2 = <P3 = ••• = 92m = (P 
W = (0,0,...,0) =0'. 



(7.31) 



As before, in the next run of the experiment we choose the same phases but we change 
the role of observers such that in the second run the first one chooses 0' while the rest of 
them choose 0, in the third run the third one chooses 0' while the rest of them choose 0, 
etc. Again, the value of the correlation function for each experiment is the same 
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2m +1 



exp 



qmW , • • • , 
2m7T 
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2m + 1 
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+ exp 



2m7T 



2m + 1 
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2m + 1 



2m exp —i 



+ exp t 



2m-K 
2m + 1 



4m^7r 
2m + 1 



exp I 



4m^7r 
2m + 1 



7^" 



(7.32) 



Using ( [7.32|) , the structure of the hidden variables correlation function ( [7.15|) and multiplying 
the above equations by each other we arrive at 



2m+l 



2m+l 



2m+l 



n 4(0', A) = n 4(0, A)- 



2m 



n 4(0, A) 



(7.33) 



k=l 



k=l 



k=l 



for every A. Because Eqm{4>' , . . . , 0') = 1 one must have 



which gives 



2m+l 

n 4(0', A) 

k=l 



2m+l 

n 4(0, A) = 1 

k=l 



(7.34) 



for every A. Thus, local hidden variables imply the following perfect correlation 

Eqm{4>, ...,0) = 1, 



(7.35) 



(7.36) 
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which is untrue because quantum mechanics gives 

Eqm{<P, 

2m + 1 L ^ 



2m + 1 

I 71 



2m + 1 



/ (2m + l)2m ^ 

+ exp i — ; 71 

\ 2m + 1 ; 



-2m + 1 
2m + 1 ' 



-2m+l 
2m+l ' 



(7.37) 



which in the hmit of 



Therefore, for each m one obtains the untrue identity 1 = 
m — > oo becomes 1 = —1. 

Now, let us see what happens when one has the even number of quA^its and the multiports 
A'" = 2m (m > 2). As before we must find appropriate phases. Let us make the following 
choice: 



= (0, 
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2m — 1 ' 2m — 1 ' 



(2m - l)7r, 
2m- 1 ' 



= 



One easily finds that 



Eqm{(, 
_ 1 

~ 2m 



+ exp i 



^2 = (P3 = ■■■ = VN-1 

^Ar = (O,O,...,O) = 0' 



, . . . ,0') = Eqm{4>', 

f 2771-1 

(2m — 1) exp —7 TT 

^ ^ ^ V 2m - 1 

= -1 



(7.38) 



Eqm{(P, . . . ,0',0) 



(2m 



2m- 1 



-7T 



(7.39) 



and 



Eqm{(I>, -,0) 
1 



2m 
1 
2m 



(2m — 1) exp 
(2m — 1) exp 



2m 



-7- 



2m - 1 
2m 



71 + exp 7 



2m(2m + 1) 
2m + 1 



-TT 



2m - 1 

Applying the same reasoning as above one has 



TT +1 
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2m 



(2m — 1) exp ( —i 



2m 



2m - 1 



TT +1 



(7.40) 

(7.41) 
oo becomes 



Again, for each m > 2 one has a contradiction, which in the limit of m - 
1 = -1. 

According to the old quantum rules of thumb one approaches the classical limit with the 
growing quantum numbers. The above results (for both the odd and the even case) once 
more show the contrary behaviour. 
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7.6 Conclusions 



The derived paradoxes can be divided into two groups. To the first group belong the series of 
paradoxes where one has more observers than input ports in the unbiased multiport beam- 
splitters. Within this group the contradiction manifests itself in different perfect correlations 
predicted by quantum mechanics and local realism. As a special case one obtains the original 
GHZ paradox. 

In the second group one has the series of paradoxes in which the number of observers 
and input ports in the unbiased multiport beamsplitters is the same. Within this group the 
derivation of the paradoxes relies on the perfect correlations (as in the first group) but the 
final result is different. The resulting contradiction between quantum mechanics and local 
realism manifests itself in the fact that local realism predicts a certain perfect correlation 
whereas quantum mechanics does not. 

An interesting feature of paradoxes within the second group is that they naturally split 
into two parts: the even and the odd number of observers (input ports). In each part the 
final contradiction has different numerical values. Furthermore, for the case of two observers 
one cannot derive the GHZ paradox with the method presented here. 

All paradoxes do not vanish with the growing dimension N of the Hilbert space describing 
each quiVit. 

In conclusion we state that the multiport beam splitters, and the idea of value assignment 
based Bell numbers, lead to a strikingly straightforward generalisation of the GHZ paradox 
for entangled quA^its. These properties may possibly find an application in future quantum 
information and communication schemes (especially as GHZ states are now observable in 
the lab ||Bouwmeester99|| ) . 



Part III 

Extension of Bell Theorem 
numerical approach 



In the present part of the dissertation a novel approach to the Bell theorem, via numerical 
linear optimisation, will be presented. This approach enables one to first of all solve the old 
question^ of generalisation of the Bell theorem to 

• more possible settings at each side of the experiment (for two and three qubits) 

• entangled pairs of quA^its (A^ > 2) (in earlier literature called spins higher than i) 

The early Bell theorem involved only two entangled qubits and two observers, who could 
choose among two local mutually incompatible observables. Generalisations involving more 
than two observers were considered by many authors and here in the previous sections. 
However, one can ask the following question what the necessary and sufficient conditions for 
local realistic description of situations in which each observer can choose from more than 
two observables, and for the case when two observables are not of a dichotomic nature are. 

There were several attempts to unite Bell inequalities for such problems, however only the 
stage of necessary conditions for local realism was reached. No one has been able to construct 
the full set of Bell inequalities for such problems which would provide the sufficient condition 
for local realistic description. 

This was a serious issue, since necessary conditions are per se weaker than sufficient ones. 
Several papers were written discussing the exploding algebraic difficulty in finding sufficient 
sets of Bell inequalities ||Peres99|] . For example Peres gives as a number of the so called 



Farkas vectors (effectively coefficients in generalised Clauser-Horne inequalities) for a two 
quiVit experiment with only two observables (non-degenerate) allowed to be measured at 
each side |[Pereslj9|| . It reads Nf{N) = {2^ -2)\ i.e., one has Nf{2) = 16, A^f(3) = 1296 etc. 



This causes that the computational time grows extremely fast not allowing one to calculate 
anything in reasonable time. 

However if one asks directly the question whether the probabilities given by quantum 
mechanics are describable via a local realistic theory, it turns out that this question can 
be formulated as a typical linear optimisation problem. Since excellent numerical methods 
exist it turns out that for the number of settings of the order of 10 for two qubits and of the 
order of 4 for three qubits, and for quA^its with up to 10, such problems are solvable on 
a standard work-station with the computing time of the order of one day. The algorithm to 
tackle these problems together with the results will be presented here. 

The threshold noise, or in other words visibility, below whichQ there is no local realistic 
model for two and three maximally entangled qubits is of the same value as in the case 



See for instance | Mermin80| , [Ardehali9l| , |W6dkiewicz91| , |Gisin92| , W6dkiewicz94 , Gisin9£, Larsson99 



In terms of the visibility above which. 



of standard Bell inequalities (involving only two settings at each side). Nevertheless, the 
presented approach gives a new method of a direct analysis of experimental data. Such data 
can be directly tested for a possibility of a local realistic model of them (without any of the 
usual hypothesis about the curves best fitting to the data). 

In the case of the second problem, entangled quA^its, new surprising results have been 
obtained. The computer data reveal growing discrepancy of the quantum predictions with 
local realism with the increasing dimensionality of the Hilbert space of the sub-systems. 
Earlier approaches suggested contrary behaviour (see , for instance, | Mermin80| , MerminS"^ ) . 



Chapter 8 

Necessary and sufficient conditions to 
violate local realism for two 
maximally entangled qubits- extension 
to more than two local settings [7] 



8.1 Introduction 



It is a common wisdom in the quantum optical community that the threshold visibility of 
the sinusoidal two-qubit interference pattern beyond which the Bell inequalities are violated 
is (for the case of perfect detectors) (see, e. g. ||Clauser78[| ). Most of the experiments 
exceed that limit (with the usual "fair sampling assumption" ) [Preedman72| , Pu88| , [Rarity9d| , 
Kwiat9^j| , [Fittman^ , [l'ittel98|] . Some difficulties to reach this threshold were observed in the 
very early experiments fClauserT^] , as well as in some recent ones involving novel techniques. 



Thus far, in atomic interferometry EPR experiments ||Hagley97|] and for the phenomenon of 
entanglement swapping ||Zukowski95| , Pan98|| , the resulting visibility is less than the magic 



71%. 

It is also well known that the Clauser-Horne inequality and the CHSH inequality are 
not only necessary conditions for the existence of local realistic models but also sufficient 
b''ine82[| (in the case of the CHSH inequality, this requires some simplifying assumptions 



3arg87|| ). However, the sufficiency proofs used involve only two pairs of settings of the local 



macroscopic parameters (e.g. orientations of the polarisers) that define the measured local 
observables. Thus, the constructions are valid for precisely those settings and nothing more, 
and there is no guarantee that the models can be extended to more settings. Consequently 
one may ask what is the maximal visibility for a model applicable to all possible settings of 
the measuring apparata, that returns sinusoidal two qubit interference fringes. It is already 
known that for perfect detectors, this value cannot be higher than or lower than 2/7r 
(this is the visibility of the recent ad hoc model by Larsson |[Larsson99[| ; for earlier models 
returning visibilities of 50% see e.g. |[W6dkiewicz91|| ). 



The knowledge of the maximal visibility of sinusoidal two-qubit fringes in a Bell-type 
experiment that still can be fully modelled in a local realistic way, may help us to distinguish 
better between 'local' and 'nonlocal' density matrices. For two two-state systems one can 



find precise conditions which have to be satisfied by density matrices describing the general 
state, pure or mixed, of the full system, that enable violation of the CHSH inequalities 
H.orodecki95||. States fulfilling such conditions are often called "nonlocal". However, since 



the CHSH inequality is necessary and sufficient only for two pairs of settings, it is not 
excluded a priori that some states that satisfy such inequalities for all possible sets of two 
pairs of local dichotomic observables, nevertheless give predictions that in their entirety 
cannot be modelled by local hidden variables. Such models must first of all reproduce the 
full continuous sinusoidal variation of two-qubit interference fringes, as well as the other 
predictions. 

It is clear that the full solution of the question would require a construction, or a proof 
of existence, of local hidden variable models which return sinusoidal fringes of the maximal 
possible visibility, that are applicable for all possible settings of the measuring apparata. 
Since this seems to be very difficult, we chose a numerical method of point wise approximation 
at a finite number of settings at each side of the experiment. Due to exponential growth of 
the computation time when the number of settings increases, we managed to reach up to 10 
settings on each side, i.e. up to 100 measurements points (which due to a certain symmetry, 
about which we will say more later, effectively can be transformed into 20 x 20 = 400 
points). The exponential growth hinders any substantial increase in this number. Such 
numerical models cannot give a definite answer concerning the critical visibility of sinusoidal 
fringes, however our calculations enable us to put forward a strong conjecture that this value 
must be indeed (see below). The numerical method presented in this chapter is applied 
to the correlation function for two qubits. The alternative approach, much better suited for 
general problems is to apply the method to probabilities of all events that are observed in 
the experiment. However, for maximally entangled two qubit state, with the interference 
visibility reduced by some noise (the generic problem studied here) application of numerical 
approach to correlation function is equivalent to a similar procedure with probabilities, and 
what is important due to the fact that it leads to a less complicated linear optimisation 
problem it allows to have numerical results for much more settings on each side than the 
method involving probabilities. Nevertheless the algorithm involving probabilities was also 
written, and it returns identical results. It will be presented in the chapter on two entangled 
quA^it problem. 

Experimentally our problem can be formulated in the following way: the two qubit state 
produced by the source does not allow for single qubit interference, and in the experiment 
less-than-perfect two-qubit fringes are obtained, due to some fundamental limitations (like 
those present in the case of entanglement swapping, e.g. [[Zukowski95 , Pan98 |) or due to 



imperfections of the devices. What is the critical two qubit interference visibility beyond 
which the observed process falsifies local realism? We shall ask these questions assuming, for 
simplicity, perfect detection efficiencies, which is possible theoretically, and experimentally 
thus far amounts to the usual "fair sampling assumption" . 

Furthermore, the problem may be investigated without the use of the assumption that the 
observed fringes are of a sinusoidal nature even though the observed two-qubit fringes in ex- 
periments with high photon counts follow almost exactly the sinusoidal curves |[Freedman72| , 
Ou88| , [RarityQO] , P:lwiat93l , [Pittman95| , [l'ittel98|| . In experiments with lower count rates, still 



with relatively good level of confidence the recorded data have approximately the same char- 
acter, and it is customary to fit them with sinusoidal curves. It is now a standard procedure 
to perform the two-qubit interference experiments by recording many points of the inter- 
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ference pattern, rather than stabihsing the devices at measurement settings appropriate for 
the best violation of some Bell inequality. Further, in some of the experiments, e.g. those 
involving optical fibre interferometers it was not possible to stabilise the phase differences 
and what is observed is just the interference pattern changing in time, and the visibility of 
the sinusoidal two-qubit fringes is used as the critical parameter |[Tapster9"^] . 

Even though the numerical calculations presented here only reach 20 x 20 points, this is 
more than enough in comparison to the experimental data. The usual experimental scans 
rarely involve more than 20 points. Further, our algorithm can be applied directly to the 
measurement data, and in that way one can even avoid the standby hypothesis that the 
fringes follow a sinusoidal pattern. The algorithm can directly answer the question: are the 
data compatible with local realism or not? Since physics is an experimental science, the 
questions about Nature get their final answers solely in this way. 



8.2 Quantum mechanical and local realistic description 
of the gedanken experiment 

We consider the following gedanken experiment. Two observers a and b perform measure- 
ments of observables a ■ a (observer a) and b ■ a (observer b) on the state defined by the 
equation (|2.1|). Observer a can choose between Na settings of the measuring apparatus, 
which are defined by vectors a, (i = 1, . . . , No) whereas observer b can choose between iVf, 
settings of the measuring apparatus defined by vectors bj (j = 1, . . . , Nh). Therefore, they 
perform Na x Nf, mutually exclusive experiments. For each experiment one has the following 
quantum probabilities (see the equation ( |6.1|) ) 

Pqlii^'nT'laify = ^(1 - /mVaai ■ bj) (8.1) 

and the quantum correlation function 

E^\,{dJ,) = -V2d,-bj. (8.2) 

There is no one qubit interference, i.e., PQ%[{l\ai) = Pqliimlbj) = \ for all i, j. 

If one assumes that the unit vectors which define the measured observables are always 
coplanar the correlation function can be simplified to i?Q^(ai, (3j) = —V2 cos (oj + (with 
the obvious definition of a, and Pj). 

Please note that in the considered here gedanken experiment there is one-to-one equiva- 
lenceQ between the quantum correlation function and quantum probabilities (see ( p.l| ) and 
( ^.2|) ). Therefore, we can use either the quantum correlation function or the quantum prob- 
abilities to describe the experiment. However, the description of the experiment in terms 
of the quantum correlation function is more convenient from the numerical point of view, 
which we will be seen further. 

^ This is possible due to "well" defined values which we ascribe to the results of measurements (here 
±1) and symmetries exhibited by the quantum probabilities. In the case of quA'^its with = 3 we will see 
that the Bell number assignment (see the chapter on GHZ paradoxes for quA^its) also enables us to use the 
quantum correlation function (|7.19D instead of the probabilities. 
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From the numerical values of the quantum correlation function we form a certain matrix 
Qy^ of quantum predictions with the entries: Qij{y2) = E^^^ifli^ bj). 

Within the local hidden variables formalism the correlation function must have the fol- 
lowing structure 

ELHv{a^,bJ) = J dXp{X)A{a,,X)B{bj,X), (8.3) 
where for dichotomic measurements 



A{di,X) = ±l (8.4) 

and 

B(bj, X) = ±1, (8.5) 

and they represent the values of local measurements predetermined by the local hidden 
variables, denoted by A, for the specified local settings. This expression is an average over 
a certain local hidden variables distribution p(A) of certain factorisable matrices, namely 
those with elements given by Mjj(A) = A{ai, X)B{bj, A). The symbol A may hide very many 
parameters. However, since the only possible values of A{ai, X) and B{bj,X) are ±1 there 
are only 2^" different sequences of the values of (A(ai, A), A{d2, A), Aidp^^, A)), and only 
2^'' different sequences of (-B(&i, A), -8(62, A), B{b]^^^, A)) and consequently they form only 
2Na+N^ matrices M(A). 

— * 

Therefore the structure of local hidden variables models of EiuviO'i-, bj) reduces to discrete 
probabilistic models involving the average of all the 2^""'"^'' matrices M(A). In other words, 
the local hidden variables can be replaced, without any loss of generality^, by a certain pair 
of variables k and / that have integer values respectively from 1 to 2^° and from 1 to 2^\ 
To each k we ascribe one possible sequence of the possible values of A{di, A), denoted from 
now on by A{di, k), similarly we replace B{bj, A) by B{bj, I). With this notation the possible 
local hidden variables models of the correlation function EiHy{di,bj) acquire the following 
simple form 

2Na 2^b 

ELHviaJ,) = J2Y.PklMS^,k)B{b„l), (8.6) 

k=l 1=1 

with, of course, the probabilities satisfying pki > and Y.t=i J2'i=i Pki = ^■ 

The special case that we study here enables us to simplify the description further. To 
satisfy the additional requirement that the local hidden variables model returns the quan- 
tum prediction of equal probability of the results at the local observation stations, that 
PQj^{l\d) = PQlj{m\b) = |, one can use the following observation. For each fc, there must 
exist a, k' ^ k with the property that A{di, k') = —A{di, k), and similarly for each /, there 
must exist an /' 7^ Z for which B{bj, I') = —B{bj, I). Then A{di, k)B{bj, I) = A{di, k')B{bj, I'), 
and thus they give exactly the same matrix of local hidden variables predictions. By assum- 
ing pki = Pk'i' the property of total randomness of local results will always be reproduced 
by the local hidden variables models, and the generality will not be reduced since the con- 



tributions of pki and pk'i' to (3^) cannot be distinguished. Hence, we will always take only 



one representative of the two pairs, reducing in this way the number of probabilities and 



See, for instance, |Wigner7C, Belinfante73 
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matrices of LHV predictions in (|8.6| ) by a factor of two^ This is equivalent to taking only 
half of different B{bj,l) (or A{ai,k) for the situation is symmetrical), say, first 2^*'"^ ones. 
Thus, ( ^.6|) acquires the form 



2Na 2^^6-1 

ELHv{ai,bj) = Y, PkiMai,k)B(bj,l), (8.7) 

k=l 1=1 

where we still assume that J2kl^i ^"1=1 ^ Pm = ^■ 

Another reduction by a factor of four is given by the fact that in the coplanar case, the 
choice of the settings may be limited on each side to ranges not greater than vr (i.e. (p < ai < 
(p + Tiy'ip < Pj <■?/' + vr). This is due to the simple observation that a model of the type ( WM 
once established for such settings can be easily extended to settings = aj + vr, = Pj + n 
by putting A{a[j) = —A{aiJ) and B{PjJ) = —B{PjJ). 

The conditions for local hidden variables to reproduce the quantum prediction with a 
final visibility V2 can be simplified to the problem of maximising a parameter V2 for which 
exists a set of 2^""'"^''"^ probabilities pki, such that 

2Na 2^6-1 

Y PkiA{a„ k)B{b„ I) = Q,,{V2). (8.8) 

k=l 1=1 

Because, for the given local settings ( ^.8|) imposes linear constraints on the probabilities and 
the visibility^, and we are looking for the maximal V2, the problem can be solved by means 
of linear programming- a certain method of optimisation. 



8.3 Linear programming and Downhill Simplex Me- 
thod 

Let us briefly describe the idea of linear programming sending more interested readers to the 
excellent book by Gass ||Gass75|| . 



The set of linear equations ( ^.8|) constitute a certain region in a D = 2^""^^''"^ + 1 dimen- 
sional real space- 2^'^^^''^'^ probabilities plus the visibility. The border of the region consists 
of hyper planes each defined by one of the equations belonging to ( p.8|) ; thus, if the equations 
do not contradict each other the region is a convex set with a certain number of vertices. 
On this convex set we define a linear function (cost function) /(pi, ■ ■ ■ ,p2Na+N^-i, V2) = V2, 
which maximum we seek. 

The fundamental theorem of linear programming states that the cost function reaches its 
maximum at one of the vertices. Hence, it suffices to find numerical values of the cost function 
calculated at the vertices and then pick up the largest one. Of course, the algorithmic 
implementation of this simple idea is not so easy for we must have a method of finding the 
vertices for which the value of the function continually increases^] so that the program reaches 

■^The method involving directly probabilities of all pairs of events avoids this artificiality. However, the 
price paid is a much longer computational time. Nevertheless, we have performed also such a calculation 
(see the last chapter). Of course the results of the presented calculation were confirmed. 

■*The visibility should also fulfil V2 < 1. 

^We must also know how to find a starting vertex. 
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the optimal solution in the least possible number of steps. Calculating the value of the cost 
function at every vertex would take too much time as there may be too many of them. 

There are lots of excellent algorithms which solve the above optimisation problem. Here 
we have used the algorithm invented by Gondzio | |Gondzio95[| and implemented in the com- 



mercial code HOPDM 2.30 (Higher Order Primal-Dual Method) written in C programming 
language. 

However, finding the maximal visibility for the given local settings of the measuring 
apparata is not enough. We should remember that our main goal is to find such local setting 
for which the threshold visibility is the lowest one. The maximal visibility V^"'^ returned 
by the HOPDM 2.30 procedure depends on the local settings entering right hand side of 
( ^.8| ). Thus, returned V.J^"'^ can be treated as the many variable function, which depends 
on Na + Nh angles in the coplanar case and two times more in the non coplanar one, i.e.. 

Hence, we should also have a numerical procedure which finds the minimum of V^"'^. 
Because we do not know much about the structure of V-J^""^ as a function of the local settings 
the only reasonable method of finding the V^""^ minimum is the Downhill Simplex Method 
(DSM) |[Nelder65[| . The way it looks for the extremum is the following. If the dimension of 



the domain of a function is Dim the DSM randomly generates Dim + 1 points. This way it 
creates a starting simplex, which vertices are these points. Then it calculates the value of 
a function at the vertices and starts exploring the space by stretching and contracting the 
simplex. In every step when it finds a vertex where the value of the function is lower then in 
others it "goes" in this direction. For more elaborate discussion of this issue see |Nelder65 



8.3.1 Numerical difficulties 

Obviously there are some numerical limitations of our computer program. First of all it is 
obvious that the time needed for computation grows with the number of local settings N^, Nh. 
For the given Na, N^ we have 2'^'""'"^''"^ -|- 1 unknowns and Na x Ni, + 2 linear constraints 
imposed on them. Furthermore, the dimension of the domain of V.^""^ is in general 2Na + 2Nh, 
which makes it harder for the DSM procedure to find a global minimum of Vg™"^. 

Secondly, the amount of the memory needed to store "hidden" matrices grows exponen- 
tially with Na and A*";,. All these problems taken together has limited our research to the 
cases where Na + Nt, < 20. 

It is worth mentioning at this point about the advantages of using the quantum correlation 
function instead of quantum probabilities as the description of the gedanken experiment. If 
we are to use probabilities^ we will have 2^"+^'' + 1 unknowns (2-^'"+^'' probabilities plus the 
visibility V2) and 4 x Na x Nb + 1 constraints. Because we struggle with the memory capacity 
and the computation time every reduction in the number of unknowns and constraints is a 
great advantage. 

Recently A. Peres [[Peresyy]] has discussed the algorithms which search for so-called Farkas 



vectors, which in turn define coefficients in generalised Bell-inequalities, the set of which is 
a sufficient and necessary condition for classical probabilistic model (here, essentially, local 
realistic) to reproduce a certain set of probabilities for pairs of experiments. However, his 
method explodes numerically much much faster then ours. Simply, our method is applied 

^For the details concerning this approach see the chapter about pairs of entangled quA^its. 
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directly to a certain finite set of specified quantum probabilities^. Whereas inequalities based 
on the Farkas lemma apply to all possible sets of probabilities. 



8.4 Results 

First of all, we have performed calculations in the coplanar case. We have checked situations 
in which the number of settings for observer a and b is the same, i.e., A^a = ^6 = 2, 3, 4, . . . , 10 
as well as the cases in which the number of settings for observer a and b is different, i.e., 
Na 7^ Nb. Because of mentioned limitations of the computer program we could only reach 
the limit of Na + N^ = 20 in both cases, i.e., for the equal and unequal number of the settings 
at each side. 

Because of the richness of Na + Nb dimensional space (coplanar case) to find a global 
minimum of V^""^ for the cases where Na + Nb < 12 we have run the DSM procedure 30 times 
with varied starting points. For the higher dimensional ones, i.e., where 13 < Na + A*";, < 15 
we have run the DSM only a few times whereas for the extremal cases, i.e., where 16 < 
Na + Nb < 20 we have not used the DSM but we simply have calculated V{^"''^ on some sets 
of angles containing (or not) the Bell angles (the definition of the Bell angles is given below). 

For Na+Nb < 15 the DSM has always converged to V^""^ = ^ within the given numerical 
precision of computation. Furthermore, within the angles Oj, Pj {i = 1 . . . Na,j = 1 . . . A^^) 
for which the optimum has been achieved there has been always the subset of four ones (we 
call the angles belonging to this subset the Bell angles) for which the maximal violation of 
the CHSH inequality occurs. In other words it means that in the quantum matrix Q (matrix 
of results for the visibility equal one) with entries calculated for the optimal local settings 
found by the DSM, i.e., settings for which V^"'^ reaches its minimum, a 2 x 2 sub matrix 13 
appears with moduli of all its elements equal to 4^, and three of them of the same sign, e.g. 



1 


1 


V2 




1 




V2 





(8.9) 



Of course, the elements of B can be scattered throughout the matrix Q but by relabelling 
the indices it is always possible, for instance, to have B in the left hand corner of Q. 

For the combinations of the number of local settings at each side of the gedanken exper- 
iment such that IQ < Na + Nb < 20 we have calculated V-J^""^ on the sets of local settings 
including the Bell angles with exactly the same result, i.e., V^™"^ = -1= whereas for the sets 

of the settings not including the Bell angles V.^^""^ > 

We have also tested non coplanar settings with exactly the same result. The visibility 
was higher than ^ if among the local settings were no the Bell angles (no sub matrix B) 

and it was exactly ^ otherwise (sub matrix B present). 
8.4.1 Exemplary numerical model 

Let us show an example of the computer solution for the 3x3 coplanar case with the 
optimisation over local settings of the measuring apparata (the DSM has been used). The 

^ Instead of the probabilities we use the correlation function but this does not change anything because 
of the mentioned equivalence of these two ways of description of the considered gedanken experiment. 
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threshold visibihty for the optimal settings found by the program is V^" 
the quantum matrix is given below 



C^-ymax 



■^max 



Q 



0.107 
-0.707 
-0.707 



-0.994 
-0.707 
0.707 



-0.399 
-0.964 
-0.266 



0.707 whereas 



^.10) 



The computer hidden variable model reproducing the above matrix reads 



-1-1 
p(+l + 1-1 
p(+l + 1-1 
p(-l-l + l 
p(+l -1 + 1 
p(+l -1 + 1 
p(-l + l-l 
p(+l -1-1 
p(+l -1-1 



+ 1 
+1 
+1 
+ 1 



i + i: 
1 + 1) 
i-i: 
1 + 11 



+1 + 1 + 1) 
+1 + 1-1) 
+1-1-1) 
+1-1 + 1) 
+1-1-1) 



0.250 
0.060 
0.135 
0.157 
0.026 
0.067 
0.055 
0.217 
0.033, 



(8.11) 



where, for instance, p(+l + 1 - 


- 1;+1 - 


1 + 1) 


denotes 


factorisable matrix M 










+1 -1 






M = 


+1 -1 








-1 +1 








(8.12) 



(T denotes ordinary transposition). 

This example clearly exhibits the characteristic trait which has been mentioned already, 
namely that in the matrix Q there is a sub matrix B corresponding to the matrix obtained 
for the Bell angles- here the elements Q21, Q31, Q32, Q22- 

The numerical precision of the entries of the matrix Q and the visibility Vp""^ is in this 
example 10"^. Of course the precision can be increased when necessary but this makes the 
time of computation longer. However, in some cases we have run the program with the 
numerical precision of order 10~® with the same result, i.e., the V^""^ = 0.707107+ 10~^. 



8.5 Application to experimental data 

To apply our computer program for analysis of experimental data we replace Qij(y2) in ( ^.8| ) 
by the measured values E^^\exp), and perform the same task. If the critical V2 returned 
by the program^ is less than 1, the data cannot be reproduced by any local hidden variable 
model. Note that one even does not have to know what the settings are. 

*In this case V2 does not have the direct interpretation of visibility. Its value tells us by what factor the 
observed values of the correlation function have to be reduced, so that a local hidden variables model exists. 



7/1 



In a recent Bell-type experiment Weinfurter and Michler^ have obtained the following 
matrix of results 



Q 



exp 



-0.894 


-0.061 


0.761 


-0.851 


0.343 


0.765 


-0.625 


0.688 


0.516 


-0.251 


0.860 


0.103 


0.226 


0.921 


-0.389 


0.530 


0.651 


-0.648 


0.855 


0.323 


-0.832 


0.852 


-0.092 


-0.843 


0.785 


-0.539 


-0.638 


0.397 


-0.795 


-0.253 



(8.13) 



V 

The program gives the verdict that the values of all entries to the matrix of results have to 
be reduced by the factor of 0.796 to be describable by local hidden variables. 

Some explanation is needed here. The entries of ( |8.13| ) give the values of the correlation 
function - there were three different setting on side A of the experiment and 27 settings on 
side B, however only 10 of them are shown here. In the actual experiment only data from 
a pair of detectors were collected. To obtain the matrix we have used the usual assumption 
that E{a,(3) = 4P(+,+;a,/5) — 1. We have also renormalised the numbers of photon pairs 
counted, so that the average of the counts over approximately two periods of settings at side 
B represents the probability of | (in concurrence with the quantum prediction). To this end 
we have used data for all 27 settings on side B. 

In the recent long-distance EPR-Bell experiment the following set of values of the corre- 
lation function was obtained^: 

/ 



Q 



exp 



V 



0.960 


-0.102 \ 


0.903 


-0.375 


0.733 


-0.660 


0.479 


-0.809 


0.191 


-0.903 


-0.120 


-0.923 


-0.429 


-0.807 


-0.666 


-0.656 


-0.842 


-0.395 


-0.951 


-0.152 


-0.953 


0.171 / 



(8.14) 



This matrix has to be reduced by the factor of 0.737 to have a local realistic description. 



8.6 Conclusions 



In his paper Peres |Peres99] has conjectured that if the CH inequalities are satisfied for 



all possible subsets of two settings on one side and two settings on the other (out of all 



^The experimental results were reported in a different context in | MichlerOOj . 

-'^"The numerical values of the entries to the matrix were provided by G. Weihs (private communication). 



Na X Nb settings) then this is sufficient for a local realistic model for all Na x A*";, settings. 
Our numerical calculations strongly support this conjecture. 

Furthermore, for the situations characterised by the condition Na+Nh < 12 we have found 
the ultimate value of the visibility still admitting the existence of local hidden variables. 
Surely, we must be aware of the fact that the DSM may have not found a global minimum 
of Vg™""^', in which the case the threshold V2 would be lower than However, we think that 
this is highly improbable, which is strongly supported by calculations for Na + Nh > 12. 

As it has been shown on two examples our computer program can be also used to analyse 
directly the raw experimental data. Here the analysis in terms of the correlation function has 
been performed. However, this is not necessary. In the chapter concerning two maximally 
entangled quiVits (including also qubits) we will present a computer program entirely based 
on the quantum probabilities rather then on the correlation function. This is important 
because the one-to-one equivalence between the correlation function and probabilities is 
valid only if one assumes that certain symmetries concerning probabilities are present (see 
( ^.1|) ). In the real experiment one may expect that these symmetries are not always fulfilled. 
The program based on probabilities utterly solves this problem allowing us to analyse real 
experiments basing solely on the observed clicks of detectors. 

To conclude, the performed calculations enable us to put forward the following conjecture: 
sinusoidal two-qubit fringes of visibility up to are describable by local realistic theories. 
At this stage we are not able to give an analytic proof of the above. However, for finite 
sets of measurement points the results of data analysis with the use of our program fully 
concur with this hypothesis. One can even say that this is more than enough, since one 
cannot experimentally test exact sinusoidal nature of the two qubit fringes. The real output 
of an experiment are count rate sequences at finite number of the measurement points which 
follow sinusoidal-like pattern. To such data our program can be directly applied giving 
verdict whether the data admit a local realistic model or not. 



Chapter 9 



Necessary and sufficient conditions to 
violate local realism for three 
maximally entangled qubits- extension 
to more than two local settings 

Exactly the same question on the threshold visibility allowing local and realistic description 
for all positions of measuring apparata can be posed in the case of correlations involving three 
qubits instead of two ones. Up to now, the threshold visibility for three qubit correlations 
(GHZ correlations) is known to be V2(3) = | (in the brackets we show that we deal with the 
visibility for three qubits to distinguish it from the visibility for two qubits V2.) |[Mermin90b |. 
However, this limit has been established for the case when there are two local settings of the 
measuring apparatus at each side of the experiment and it has been obtained in a standard 
way, i.e., as a condition for the violation of an appropriate Bell inequality. 

From the previous section we know that there is a more direct method of finding the 
numerical value of the threshold visibility, which in addition gives necessary and sufficient 
conditions for the existence of local hidden variables for the given set of local settings for 
each observer. In this section we show the application of the presented numerical method 
to the GHZ correlations, i.e., to the maximally entangled state of three qubits. Surprisingly, 
the obtained results bear great resemblance to those for two qubits. 

9.1 Description of the method 

To this end let us consider the following maximally entangled state of 3 qubits 

IV^) = ^(|0)l|0)2|0)3+|l)l|l)2|l)3) (9.1) 

where is the i-th state of the j-th qubit. 

Each observer measure the observable n ■ a, where n = a,b,c {a for the first observer, b 
for the second one and c for the third one), n is a unit vector characterising the observable 
which is measured by observer n and a is a vector the components of which are standard 
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Pauli matrices. As in the case of two qubits the family of observables n ■ a covers all possible 
dichotomic observables for a qubit system, endowed with a spectrum consisting of ±1. 

The probability of obtaining the result m = ±1 for the observer a, when measuring the 
observable characterised by the vector a, the result / = ±1 for the observer 6, when measuring 
the observable characterised by the vector h and the result k = ±1 for the observer c, when 
measuring the observable characterised by the vector c is equal to 

PqAiim, /, A;; a, 6, c) = |(1 + mlaabs + mkasc^ + Ikh-iC^ 

+mlkY.l,p,s=iMrpsarbpCs), (9.2) 

where components of vectors a, 6, c and where nonzero elements of the tensor 

Mrps are Mm = 1, M122 = —1,^212 = — l,M22i = —1- In spherical coordinates vectors 
a, b, c read 

n = (cos0„ sin6'„, sin0„ sin6'„, cos6'„), (9.3) 

where On G [0, vr], G [0, 2tt]. From now on we will be considering only the measurement of 
the observables characterised by vectors with the zero third component, which is equivalent 
to putting 6n = 7t/2. Thus, the formula (|9.2|) acquires simpler form (we have replaced 
(pa, (pb, <Pc by a, /3, 7 respectively) 

PQMim,l,k;a,P,-f) = 1(1 + mlkJ2r,p,s=iMrpsarbpCs) (9.4) 

in which only a term responsible for three qubit correlations is present. 

The probabilities of obtaining one of the results in the local stations reveal no dependence 
on the local parameters, PQM(^|a) = PQM{fn\l3) = Pqm(^|7) = |- Similarly, the probabil- 
ities describing two qubit correlations do not reveal dependence on the local parameters 
either, Pqm(/, m|a, /?) = PQM{m,n\P,-f) = PQMil,n\a,-f) = \. 

As usual, if V2(3) < 1 we replace (|9.4| ) by 

P^f\mJ,k\a,P,^) = lil + mlkV2{3)Er,p,s=iMrpsarbpCs). (9.5) 

The quantum prediction for the correlation function with reduced visibility reads: 

EQM\a,f3,l) = Eln,i,k=~i^^kP{m,l,k;a,(3,^) 
= ^2(3) E'p,.=i MrpsarbpCs = V2{3) cos(a + /5 + 7). (9.6) 

and there is no single and two qubit interference^. 

Now, we proceed analogously to the case of two qubits. Observer a chooses between Na 
settings of the measuring apparata ai, . . . , oat^, observer b between A^;, settings Pi, ... , (3n^ 
and , finally, observer c between Nc settings 71, . . . ,1Nc,- For each triple of local settings we 
calculate the quantum correlation function ( |9.6| ) E^j;f\ai, f3j,'yk), where i = 1, . . . ,Na,j = 

1, . . . , Nb, /c = 1, . . . , Nc. Thus we have a tensor^ Qijk{V2{S)) = E^^\ai, f3j,'^k) of quantum 
predictions. 

^We again use the quantum correlation function instead of quantum probabilities. In the case considered 
here these two approaches are equivalent- see the discussion in the previous chapter. 

^We call it a tensor for it has three indices. This has nothing in common with the transformation 
properties of this object. 
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Within the local hidden variables formalism the correlation function must have the fol- 
lowing structure 

ELHv{au Pj, Ik) = J dp{X)A{ai, A)E(/3„ X)C{^k, A), (9.7) 
where for dichotomic measurements 

A{ai, A) = ±1 
B{f3„X) = ±l 

C(7fe,A) = ±l, (9.8) 

and they represent the values of local measurements predetermined by the local hidden 
variables, denoted by A, for the specified local settings. This expression is an average over a 
certain local hidden variables distribution p(A) of certain factorisable tensors, namely those 
with elements given by Tijk{X) = A{ai, X)B{Pj, A)C(7fc, A). Since the only possible values of 
A{ai, A), B{i3j, A) and C(7fc, A) are ±1 there are only 2^" different sequences of the values 
of (A(ai,A),...,A(aAr^,A)), 2^" different sequences of {B{(3i, X), B{(3n,, X)), 2^^= different 
sequences of (C(7i, A), C(7Ar^, A)) and consequently they form only 2'^°"'"^''+-^'= tensors 
Tijk{X). 

Therefore the structure of local hidden variables models of E^Hvicii, Pj,1k) reduces to 
discrete probabilistic models involving the average of all the 2^"+^''+^= tensors Tijk{X). 
In other words, the local hidden variables can be replaced, without any loss of gener- 
ality, by a certain triple of variables /,m, n that have integer values respectively from 
1, . . . , 2^", 1, . . . , 2^'', 1, . . . , 2^". To each / we ascribe one possible sequence of the possi- 
ble values of A{ai,X), denoted from now on by A{ai,l), similarly we replace B{Pj,X) by 
B{Pj,m) and (7(7^, A) by C{'yk,n) . With this notation the possible local hidden variables 
models of the correlation function ELHv{c(i, Pj,'lk) acquire the following simple form 

2Na 2^b 2^= 

EiHvioih Pj, Ik) = ^J2Y1 PimnA{ai, l)B{pj, m)C{'yk, n), (9.9) 

1=1 m=l n=l 

with, of course, the probabilities satisfying pi^n > and J2f^l J2m=i Ylf!^iPimn = 1- 

Not all of tensors Tijk{lmn) are different. It is easy to check that only one fourth of them 
are. The situation is similar to that with two qubits and enables one to simplify the actual 
computer program. 

For the given local settings of the measuring apparatus at each side of the experiment 
we want to find the maximal V2(3) still admitting the local realistic description in the form 
( |9.9| ). Then we want to find such local settings for which this maximal V2(3) reaches its 
minimum. From the previous section we know how to cope with such a problem. We use 
exactly the same numerical approach, i.e., the HOPDM 2.30 and the DSM procedures. 



9.2 Results 

We have checked three experimental situations: Na = = Nc = 2,3,4,5 with the result 
that the threshold visibility admitting local hidden variables is V2(3) = |. This result is in 
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concurrence with the threshold visibihty obtained earher in |[Mermin90b|| with the usage of 
appropriate Bell inequalities. 

Again, because of the complexity of the space being the domain of the V^"'^{3) function 
(defined in analogy to that for two qubits)- to find a global minimum- for the case Na = 
Ni) = Nc = 2,3 we have run the amoeba procedure 30 times with varied starting points. For 
Na = Nb = Nc = 4: we calculated ¥^""^{3) on 9000 randomly chosen sets of the local settings 
whereas in the case Na = Ni, = Nc = 5 we have calculated ¥^^""^{3) on the following set of 
the local settings: ai = 0,^2 = 7r/8,a3 = 7r/4,a4 = 37r/8,a5 = 7r/2,/5i = 7i = — 7r/4,/32 = 
72 = -vr/8,/33 = 73 = 0,/54 = 74 = n/S^/S^ = 75 = n/A. 

An interesting feature of the results is that, as for two qubits, the threshold visibility 
ymax^r^^ = | is always achieved for such settings of the measuring apparata which include as 
a subset the settings giving maximal violation of the inequalities presented in [[Mermin90b|| . 
This is in analogy with the two qubit case in which the threshold visibility of is always 
obtained if among the settings one has a subset which lead to maximal violation of the CHSH 
inequality (for the maximally entangled state). 



9.3 Conclusions 

The presented numerical approach to the three qubit GHZ correlations gives the sufficient 
and necessary conditions for the existence of local hidden variables for the given experimental 
situation, i.e., for the fixed number of positions of the measuring apparata at each side of 
the experiment. 

For the cases where Na = Nb = Nc = 2,3 we have found such numerical values of the 
local settings for which the visibility admitting local hidden variables has the lowest possible 
value. Up to the possibility that the DSM procedure has not succeeded in finding the global 
minimum of ^"^"^(3) the visibility V2(3) = | is the ultimate limit drawing the borderline 
between local hidden variables and quantum mechanics for these cases, i.e., for 2 and 3 
settings of the measuring apparatus at each side of the experiment. As far as I know these 
are the first results giving the necessary and sufficient conditions for violation of local realism 
in the GHZ caseQ 

For Na = Nb = Nc = 4: the visibility returned by the program for every random choice of 
local settings has been always higherQthan |. 

In the extremal case, i.e., for Na = Nb = Nc = 5 we have found the threshold value for 
local settings including as a subset settings giving maximal violation of Mermin's inequality 
M.ermin90B| with the result that V^"'^{3) = | (the DSM has not been used). 



Unfortunately, due to the computer time and memory limitations (of the same origin 
as in the two qubit case) we could not check more settings of the measuring apparatus. 
Nevertheless, we suppose that increasing the number of settings will not lead to visibility 
lower than 1^(3) = |. 

The important aspect of the presented analysis of the GHZ correlations is that, just like 
in the case of two qubits, the numerical approach can be directly apphed to measurement 
data. 



The inequalities found in | Mcrmin90b{ give only necessary condition for the existence of local realism. 



^Sometimes it was very close to ^ . Obviously, it is extremely difficult to find at random such local settings 
for which visibility equals exactly 2- 
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We should also mention that the program based on the quantum probabilities (the idea 
is presented in the chapter 10) has also been written for considered here three maximally 
entangled qubits. Some calculations have been performed but due to the time limitations we 
could not test so much cases as we did with the program based on the correlation function. 
As one expects, the results for the tested cases {Na = Ni, = Nc — 2,3 with the DSM 
optimisation) are the same, i.e., the threshold visibility equals |. 
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Chapter 10 



Entangled pairs of quA^its: the 
violation of local realism increases 
with N [8] 



10.1 Introduction 

John Bell has shown that no local realistic models can agree with all quantum mechanical 
predictions for the maximally entangled states of two qubits. After some years researchers 
started to ask questions about the Bell theorem for more complicated systems. The most 
spectacular answer came for multiple qubits in the form the GHZ theorem [ preenberger89|| : 
the conflict between local realism and quantum mechanics is much sharper than for two 
qubits, and can be shown even at the level of perfect EPR-type correlations. 

The other possible extension are entangled states of pairs of quA^its (3 < A^). First re- 
sults, in 1980-82, suggested that the conflict between local realism and quantum mechanics 
diminishes with growing N ||Mermin8q , |Mermin82| , |Garg82|| . This was felt to be in concur- 



rence with the old quantum wisdom of higher quantum numbers leading to a quasi-classical 
behaviour. However, the early research was confined to Stern-Gerlach type measurements 
performed on pairs of maximally entangled spins [[Mermin80| , [Mermin82i |Garg82|| . Since 



operation of a Stern-Gerlach device depends solely on the orientation of the quantisation axis, 
i.e. on only two parameters, devices of this kind cannot make projections into arbitrary or- 
thogonal bases of the subsystems. That is, they cannot make full use of the richness of the 
A^-dimensional Hilbert space. 

Wodkiewicz ||W6dkiewicz94 ] proposed to employ the measurement of the observable of a 
different kind (projection onto a coherent state) and the Clauser-Horne inequality. Even in 
this case in the limit of A^ — > oo the violation vanishes. 

In early 1990 's Peres and Gisin ||Peres92| , |Gisin92|| have shown, that if one considers 
certain dichotomic observables applied to maximally entangled pairs of quA^its, the violation 
of local realism, or more precisely of the CHSH inequalities, survives the limit of A^ — oo and 
is maximal there. However, for any dichotomic quantum observables the CHSH inequalities 
give violations bounded by the Tsirelson limit ||'lsirelson80||, i.e. limited by the factor of \/2- 



Therefore, the question whether the violation of local realism increases with growing A^ was 
still left open. 



so 



There are some reasons to suspect that violations of local realism should get stronger 
with increasing N. For systems described by observables which are at least three valued 
the Bell-Kochen-Specker theorem [ Cjleason57| , |Kochen6"7 , Bell66 | on non-contextual hidden 
variable theories can be applied. This means that any realistic theory of local observations 
must be inevitably contextual. In contradistinction, the original Bell theorem is formulated 
for subsystems for which such problems do not arise. 

In this chapter we show that violation of local realism indeed increases with growing 
if one uses non-degenerate observables- already introduced unbiased symmetric 2N ports. 

As a "measure" of the strength of violation of local realism in this chapter the threshold 
noise admixture (see ( p.l4|) ) still allowing a local and realistic description of the gedanken 
experiment will be used. Its link with the visibility is simple, namely = 1 — V^. 
However, as the results for objects belonging to Hilbert spaces of different dimensions will 
be compared, using the parameter Fn is more objective. 



10.2 Description of the gedanken experiment 

The general idea is the same as in the previous chapters. We analyse a Bell-type experiment 
with two quA^its flying towards two spatially separated observers A and B. We assume that 
the quA^its are prepared in the mixed state p{Fn) 

P{Fn) = FNpnoise + (1 " FN)Pmax, (10.1) 

where pnoise = ® I {I is an N x N identity matrix) and pmax is the projector \ip){'ip\ on 
the maximally entangled state 

1 ^ 

\i^) = ^r^\k)A\k)B. (10.2) 




Figure 10.1: The experiment of observer A and B with two maximally entangled quA^its. 
Each of their measuring apparata consist of a set of A^ phase shifters just in front of an 
2N port Bell multiport, and A^ photon detectors -D^, Di (perfect, in the gedanken situation 
described here) which register photons in the output ports of the device. The phase shifters 
serve the role of the devices which set the free macroscopic, classical parameters that can be 
controlled by the experimenters. The source produces a beam-entangled two quA^it state. 



Observer A can choose between the non commuting observables Ai,A2 and observer B 
also can choose between the non commuting observables Bi, B2 (each observable for observer 
A and B is characterised by some set of local parameters (knobs)) . We assume that the 
spectrum of each observable consists of points, which we enumerate by subsequent natural 
numbers k,l = 1,2, . . . , N, where index k refers to observer A and index / to observer B. 
Thus, the observers can perform 2x2 mutually exclusive global experiments. 

The quantum probability distribution for the specific pairs of results k and /, provided a 
specific pairs of local observables is chosen {Ai and Bj), will be denoted by Pqljik, l\Ai, Bj). 
In the case considered here the quantum probabilities read 



Pqm{^i ^Ai, Bj) 

= + (1 - FN)P^S!{k,l\A,B,), (10.3) 

where PQ^{k, l\Ai, Bj) is the probability for the given pair of events for the pure maximally 
entangled state (|10.2| ). Of course, the exact form of Pq^^k, l\Ai, Bj) depends on the specific 



choice of the observables Ai,Bj. According to quantum mechanics the set of AN"^ such 
probabilities is the only information available to the observers. 



10.3 Local realism and joint probability distribution 



It is well known (see, e. g. [^TneS^, | [Peres99[| ) that the hypothesis of local hidden variables 



is equivalent to the existence of a (non-negative) joint probability distribution involving all 
four observables from which it should be possible to obtain all the quantum predictions 
as marginalsQ. Let us denote this hypothetical distribution by Pnvik, m; I, n\Ai, A2, Bi, B2), 
where k and m, represent the outcome values for observer A observables (/ and n for observer 
B). In quantum mechanics one cannot even define such objects, since they involve mutually 
incompatible measurements. The local hidden variable probabilities Phv{---) are defined as 
the marginals 

PHv{k,l\Ai,Bi) = EmEnPHv{k,m;l,n), 
Pnvik, n\Ai, B2) = Em J2i Pnvik, m; I, n), 
PHvim, l\A2, Bi) = Efc En PHvik, m; /, n), 
PHvim, n\A2, B2) = Y.k J2i Pnvik, m; I, n), 

(10.4) 

where Pnvik, fn; I, n) is a short hand notation for Pnvik, m; I, n\Ai, A2, Bi, B2). The 4 x A^^ 
equations ( |10.4| ) form the full set of necessary and sufficient conditions for the existence of 
local and realistic description of the experiment, i.e., for the joint probability distribution 
Pnvik, TTt;l,n). The Bell theorem says that there are quantum predictions, which for F^r 
below a certain threshold cannot be modelled by ( |10.4| ), i.e. there exists a critical below 



which one cannot have any local realistic model with Pnvik, l\Ai, Bj) = PQ'^jik,l\Ai, Bj). 
Our goal is to find observables for the two quA^its returning the highest possible critical F^. 



^The simple proof of this statement is given in the Appendix A. 
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10.4 Linear programming 



The set of conditions (|10.4| ) with PQ^.j{k,l\Ai, Bj) replacing PHvik,l\Ai, Bj) imposes hnear 
constraints on the N"^ "hidden probabihties" Pnvik, fn; I, n) and on the parameter F/v, which 
are the nonnegative unknowns. We have more unknowns + than equations (4A^^), and 
we want to find the minimal for which the set of constraints can still be satisfied. Then 
we want to find such local settings characterising the observables that the found minimal 
Fm reaches its highest possible value|. From the previous chapters we know how to handle 
such a problem. We use our sledge hammer: HOPDM 2.30 and the DSM procedures. 



10.5 Observables 



In our numerical calculations we have used the observables defined by unbiased multiport 
beamsplitters (for the definition of such a device see chapter devoted to the GHZ paradoxes 
for quMts). 

The quantum prediction for the joint probability PQ^f{k,l\Ai,Bj) to detect a photon 

at the k-th output of the multiport A characterised by phase shifters (f)\ = {(f)\{i), ...0^(z)) 
{i = 1,2) and another one at the l-th output of the multiport B characterised by phase 
shifters 0-^ = (0^(j), ...0^(j)) (j = 1,2) for the maximally entangled state ( |10.2|) can be 
derived from (\IA\j with M = 2 and reads: 



p^^^ik,i\$\ji) 



N 



COS 



(10.5) 



where $^(ij) = 



+ </'b(j) + [^(^ + ^ - 2)]^. The counts at a single detector, of 
course, do not depend upon the local phase settings: PQM{k\Ai) = PQM{l\Bj) = 1/N for all 



10.6 Results 



The results are depicted in ( |10.2|) . We see that F/v continuously increases with growing N 
exhibiting opposite behaviour to the results obtained in earlier works ||Mermin80| , |Mermin8^ , 
Garg82| , |Peres92| , |Gisin92| , |W6dkiewicz9l . 

For instance, let us compare ( |10.2| ) with the results obtained in | |Gisin92| |, which are 
depicted in ( |10.3|) . We see that for = 2 both results are identical. However, starting from 
> 3 the results in ( |10.3D never exceed the 1 ~ This is due to the fact that in | |Gisin92 
dichotomic observables were used and the threshold Fn was obtained as a condition for 
violation of CHSH inequality. This behaviour is in agreement with the theorem proved by 
Tsirelson [ [IsirelsonSQ] , which states that the CHSH inequality for dichotomic observables 



can be violated only if F/v is lower than 1 — 

^Please, notice that here we are looking for the threshold value of the noise admixture Fn and not the 
threshold Vn as in the chapters concerning entangled qubits. Because Vn = 1 — Fn to maximal visibility 
refers the minimal noise admixture. 
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Figure 10.2: The threshold noise admixture Fn versus the dimension of the Hilbert space of 
quA^its. For = 3 basing on the computer calculations_we were able to find the analytical 
expression for the threshold noise admixture F3 
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Figure 10.3: The threshold noise admixture F/v versus the dimension of the Hilbert space of 



quA^its obtained by Gisin and Peres ||Gisin92 
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A few words of comment are needed. One may argue that because of a quite large number 
of local macroscopic parameters (the phases) defining the function to be maximised with the 
DSM procedure we could have missed the global minimum. While this argument cannot 
be ruled out in principle, we stress that in that case the ultimate violation would even be 
larger. This would only strengthen our conclusion that two entangled quA^it systems are in 
stronger conflict with local realism than two entangled qubits. 

An important question is whether unbiased multiports provide us with a family of ob- 
servables in maximal conflict with local realism. For a check of this question we have also 
calculated the threshold value of for the case where both observers apply to the incoming 
qutrit (A^ = 3) the most general unitary transformation belonging to a full SU(3) group 
(i.e. we have any trichotomic observables on each side). Again we have assumed that each 
observer chooses between two sets of local settings. However, in this case each set consists of 
8 local settings rather than the three (effectively two) in the tritter case. The result appears 
to be the same as for two fritters, which suggests that fritters (an perhaps generally unbiased 
multiports) are optimal devices to test quantum mechanics against local realism for A^ = 3 
(for all A^). 



10.7 Exemplary analytical model for = 3 

The presented above approach based on the idea of marginals and joint probability distri- 
bution compatible with them is the most general one, i.e., it does not assume any symme- 
tries within the tested quantum probabilities (marginals)Q However, this approach has one 
drawback. Namely, it involves lots of linear constraints imposed on the joint distribution 
probability. In the most general case, i.e., when observer A chooses Na observables Ai and 
observer B chooses A^^ observables Bj (each observable has A^ point spectrum) there are 
Na X Nb X A^^ linear constraints that must be fulfilled by the joint probability distribution 
(we does not count the constraint that < 1). This considerably lengthens computation 
time. 

On the other hand the quantum probabilities (marginals) that we consider here exhibit 
some symmetries. Therefore, one can ask the question if, as it has been with two and three 
maximally entangled qubits (see previous chapters), it is possible to use the properly defined 
quantum correlation function instead of quantum probabilities. 

In this section we show that the answer is yes but only forQ A^ = 3. Considering the 
quantum correlation function will benefit in calculating the threshold F3 for more local 
settings of the measuring apparata than with the program presented in this chapter. 

To this end let us consider two maximally entangled qutrits (quA^its with A^ = 3) and the 
quantum correlation function defined in ( [7.7D with M = 2 and A^ = 3. We remember that 
it has been obtained by ascribing to the local results of measurements subsequent powers 
of a = exp(^). We can prove that there is one-to-one equivalence between quantum 

^Therefore it can be applied for any marginals not necessarily quantum ones. 

'^We have also performed calculations for iV > 4 with the use of the quantum correlation function ( [7.7| ) 
with A/ = 2. The correlation function can be reproduced for growing A'' for lower and lower values of the 
noise fraction Fjv. Of course, this effect has nothing to do with the possibility of a local realistic description 
for such cases. Simply the correlation functions defined with Bell numbers for TV > 4 start to wash out the 
details of the full set of probabilities describing the experiment. 
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probabilities ( [7 .41) and the quantum correlation function ( [7.7|) . 

According to ( |7.4| ) for = 3 and for any local settings (pA,<pB quantum probabilities 
P'^'^ {kl\(l)A, (ps) can be divided into three groups 





= Pqai^) 


= Pqm{21) 


= PQAim 






= -Pqm(23) 


= PQKim 


p3 


= Pqm{IS) 


= Pqm(31) 


= Pqm{'22), 



(10.6) 

where we do not show the dependence on local settings to shorten the notation. This allows 
us to write the correlation function Eqm{(P\, (Pb) ^^e form (for convenience we put F3 = 0) 

Eqm = 3(P<Jm + «'^QM + «^3m), (10.7) 

where again we do not show the dependence of the correlation function on local settings. 
Using the identity | = -Pqa/ + Pqm + Pqm a + + = we arrive at 

^(Eqm - 1) = (a - 1)PIm - (2 + a)P^M, (10-8) 

which allows us to express all quantum probabilities through the correlation function (now 
F3 is arbitrary) 

1 2 
Pqm = g + qP^Eqm 

Pqm = 1- ^(v^/mEg, + ReE§^,,) 

Pqm = 1 + liVsimE^^j^ - ReE^Q'u)- (10-9) 

Therefore, we can apply the approach based on the quantum correlation function presented 
in the chapter concerning two and three qubits. 

For every measurement of the pair of observables Ai and Bj {i = 1, . . . , Na] j = 1, . . . , iVb) 
defined by the fritter and the phase shifters we calculate the quantum correlation function 
obtaining this way the quantum matrix of predictions Q with entries Qij{F^) = Eq'j^^cp^, (Ps)- 
Note that now the observer A has Na observables to choose from (B has Nb)- 

The hidden variable model in this case reads 

ELHviAJi) = PkiMA,k)B{$'^,l), (10.10) 

k=l 1=1 

where the functions A{(j)\, k) = and -B(0'b, V) = a" (m, n = 1, 2, 3). If it is to reproduce 
the quantum correlation function it must fulfil the following set of Na x A";, linear constraint^ 

QdF3)=T. PkiMA,k)B{$^s,l). (10.11) 

k=l 1=1 

The number of relevant pki can be reduced by a factor of 3 using a develop ment of the trick that has 
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led to the representation of the correlation function for two qubits in the form of (S.7). 



Using HOPDM 2.30 and the DSM we can find the threshold still admitting a local 
and realistic description of the gedanken experiment. The results are identical with those 
obtained by the previous method. However, since the algorithm is now much faster we 
were able to perform calculations involving up to 5 different phase settings on each side (i.e. 
allowing each observer to have a choice of up to 5 different observables to measure). Just as 
in the two qubit case the optimal stayed unchanged at the value obtained for the problem 
with Na = Nh = 2. 



10.7.1 Explicit model for extremal case 

Exemplary optimal settings for violation of local realism in the experiment with two qutrits 
and Na = Nh = 2 are 

0a' = (0, 0, 0) 
=(0.g.-g) 

fs' = {0,-l,^A- (10.12) 
For such settings the quantum matrix Q (with F3 = 0) read^ 

Q=(n' ? V (10.13) 



Q2 Q 



where Qi = - and Q2 = -|(1 + 2i). Please, note that the all entries of ([TaT3|) 

have the same modulus equal 

The hidden variables can only reproduce the matrix of correlation function given by 
(1 — F^)Q with F3 < i^~6v^ _ Por the threshold maximal F3 the explicit model is given by 





( a' 




) +Pl 






[ 






( a 












a j 



i^""^ ]+<![ ) = (i-i^3)Q, (10.14) 

where p = q = 

This model has to be understood in the following way. Consider the first term. The 
rank-one matrix 



(10.15) 



^We present here result of an algebraic re-calculation of the computer output. Perhaps, the most exciting 
aspect of such exercise is the exact value of F3 = ^^-^^ ^ 



on 



can be factorised into column and row matrices built out of powers of a, in the following 
three ways: 






(10.16) 



Therefore with probability p' = each of this factorisations is present in the model of Q^'-' 
(compare ([TaTOl) and ( |TOTl| )). 



10.8 Conclusions 

It is evident, that indeed two entangled quA^its violate local realism stronger than two 
entangled qubits, and that the violation increases with A^. It is important to stress that 
the values were obtained using four independently written codes, one of them employing a 
different linear optimisation procedure (from the NAG Library). 

It is interesting to compare our results with the limit ||Horodecki99|] for the non separabil- 



ity of the density matrices |[Werner89| of the two entangled systems. The fact that this limit, 
j^^, is always higher than ours indicates that that requirement of having local quantum 
description of the two subsystems is a much more stringent condition than our requirement 
of admitting any possible local realistic model. 

It will be interesting to consider within our approach different families of states, general- 
isations to more than two quA^its, extensions of the families of observables, to see if a wider 
choice of experiments than can be performed on one side (i.e., more than two) can lead to 
even stronger violations of local realism, and finally to see experimental realizations of such 
schemes. 

Finally we should mention that this approach can be extended to the research into the 
threshold efficiency of detectors for such experiments. The computer program for imperfect 
detectors has been already written. 



Appendix A 

Proof of equivalence of the existence 
of local hidden variables and a joint 
probability distribution for 
incompatible measurements 



This appendix is a re-derivation of known results (see e.g. |[l:''ine82|| ). 



We consider the experiment with two observers A and B. Observer A can choose to 
measure Na observables Ai {i = 1,2,..., Na) whereas observer can choose to measure Nb 
observables Bj (j = 1, 2, . . . , Nb)- In the measurement of the observable Ai observer A can 
obtain N results, which we denote by ki [i = 1,. . . ,N). Similarly, in the measurement of 
the observable Bj observer B can obtain N results, which we denote by Ij [j = 1, . . . ,N). 
In the experiment the observers can only measure joint probabilities P{ki,lj\Ai, Bj), i.e., 
probabilities of obtaining the result when measuring the observable A^ and the result Ij 
when measuring the observable Bj. 

First, let us show that the existence of stochastic local hidden variables recovering the 
marginals P{ki,lj\Ai, Bj) implies the existence of a joint probability distribution 

PHv{ki, . . . , kNyi, h, ■ ■ ■ , ^NbI^I: • • • ; ^Na: ^i, . . . , B^g) (A.l) 

compatible with the given marginals. The existence of local hidden variable space A and 
the "hidden" probabilities Pi{ki\X,Ai) and Pj{lj; \Bj) (see (|2.13|) ) allows one to define the 
joint probability distribution ( [A.l| ) in the following way 



Pnviki, . . . , /ctv^, /i, . . . , l]\fj^\Ai, . . . , An^, Bi, . . . , Bj^^] 

Na Nb 



[ dp{x) n p^{h\x, A,) n p^m, b,). (a.2) 



i=l 3=1 



It is evident that ( |A.2D returns as marginals P(A;j, 5^). 

Now let us assume that ( |A.1| ) exists and is compatible with the marginals P{ki, lj\Ai, Bj] 



We can define the local hidden variables as follows. To each sequence of possible results 
{ki, ^2, . . . , ^AT^, /i, ^2, ■ • ■ , Inb) we ascribe a hidden variable, which we denote as A(fci, ^2, . . . , 
k^A^ h, h, ■ ■ ■ , ^Nb)- Therefore, the local hidden variable space A consists of iV^-^+^s different 



no 



hidden variables. On this space we define the discrete probabihty distribution of local hidden 
variables p{X) in the following way 

p{\{ki, k2, . . ■ , kj^^, li,l2, ■ ■ ■ , ^Nb)) 

= Pnviki, . . . , /ctva! h, ■ ■ ■ , ^iVsl^i, • • • , ^Na^ Bi, . . . , B^g). (A. 3) 

The hidden probabilities Pi{ki\\,Ai) and Pj{lj\X, Bj) we define as 

Pi{ki\X{ki, k2, . . . , kj^^, li, I2, ■ ■ ■ , ^Afs), Ai) = S-^, ,^, 

Pj{lj\X{ki, h,---, kNA, h, h, ■■■Jnb)^ Bj) = 5;-^, (A.4) 

where S denotes the Kronecker delta function. Then the marginals are recovered in the 
following way 

P{ki, lj\Ai, Bj) 

= piHh, ^2, • • • , kNA, h, • • • , lNB))^k,,kAi,h^ (^-^^ 

A 

where we sum over all A^^^+^s hidden variables A. Using the above definitions ( |A.5| ) can be 
rewritten as 

P{ki,lj\Ai,Bj) =J2pi\)Piih\X,A,)P,{l,\\,Bj), (A.6) 

A 

i.e. we have a typical structure of a local hidden variable model (for hidden variables forming 
a discrete set). 
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